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Progress  06/2011-05/2012 
Overview 

In  this  project  we  aim  to  directly  realize  a  model  system  of  strongly  correlated 
electrons  moving  in  two  dimensions  using  ultracold  fermionic  atoms  stored  in  a 
sheet  of  light.  The  goal  is  to  create  high-temperature  superfluids  in  two 
dimensions,  to  establish  interferometry  and  magnetometry  with  these  systems, 
and  to  study  the  phase  diagram  of  two-dimensional  Fermi  gases  with  arbitrary 
interactions  and  spin  imbalance.  Fast  rotation  will  mimic  high  magnetic  fields  and 
allow  the  approach  to  the  Quantum  Flail  regime.  The  system  shares  traits  with 
FHigh-Tc  materials,  where  super  currents  flow  between  stacks  of  weakly  connected 
two  dimensional  planes.  This  year  has  seen  rather  spectacular  progress:  We 
could  study  the  evolution  of  fermion  pairing  from  three  to  two  dimensions; 
determine  the  equation  of  state  of  strongly  interacting  Fermi  gases  across  the 
superfluid  transition;  produce  mixtures  of  sodium  and  fermionic  potassium  for  the 
study  of  impurities  in  lower  dimensions;  and  create  and  study  spin-orbit  coupling  in 
a  degenerate  Fermi  gas. 


1 .  Evolution  of  Fermion  Pairing  from  Three  to  Two  Dimensions 

Ariel  T.  Sommer,  Lawrence  W.  Cheuk,  Mark  Jen-Flao  Ku,  Waseem  S.  Bakr,  Martin  W. 

Zwierlein 

Phys.  Rev.  Lett.  108,  045302  (2012) 

Plighlighted  as  a  Viewpoint  in  Physics  5,  10  (2012)  by  Mohit  Randeria 

Interacting  fermions  in  coupled  two-dimensional  (2D)  layers  present  unique  physical 
phenomena  and  are  central  to  the  description  of  unconventional  superconductivity  in  high- 
transition-temperature  cuprates  and  layered  organic  conductors.  Reduced  dimensionality 
enhances  the  effect  of  fluctuations,  while  interlayer  coupling  can  stabilize 
superconductivity  and  even  amplify  the  transition  temperature.  A  fermionic  superfluid 
loaded  into  a  periodic  potential  should  form  stacks  of  two-dimensional  superfluids  with 
tunable  interlayer  coupling,  a  key  ingredient  of  the  model  proposed  by  Anderson  to 
explain  high  transition  temperatures  observed  in  the  cuprates.  For  deep  potentials  in  the 
regime  of  uncoupled  2D  layers,  increasing  the  temperature  of  the  gas  is  expected  to 
destroy  superfluidity  through  the  Berezinskii-Kosterlitz-Thouless  mechanism,  while  more 
exotic  multi-plane  vortex  loop  excitations  are  predicted  for  a  3D-anisotropic  BCS 
superfluid  near  the  critical  point. 

In  this  work,  we  studied  fermion  pairing  across  the  crossover  from  3D  to  2D  in  a  periodic 
potential  of  increasing  depth.  We  follow  the  evolution  of  fermion  pairing  in  the  dimensional 
crossover  from  3D  to  2D  as  a  strongly  interacting  Fermi  gas  of  6Li  atoms  becomes 
confined  to  a  stack  of  two-dimensional  layers  formed  by  a  one-dimensional  optical  lattice. 
Decreasing  the  dimensionality  leads  to  the  opening  of  a  gap  in  radiofrequency  spectra, 


even  on  the  BCS-side  of  a  Feshbach  resonance.  With  increasing  lattice  depth,  the 
measured  binding  energy  EB  of  fermion  pairs  increases  in  surprising  agreement  with 
mean-field  theory  for  the  BEC-BCS  crossover  in  two  dimensions. 
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Figure  1  Figure  1  Evolution  of  Fermion  Pairing  from  Three  to  Two  Dimensions.  Radio-Frequency 
Spectra  show  the  opening  of  a  pairing  gap  as  the  Fermi  Gas  is  more  and  more  confined  to  two 
dimensions.  VO  denotes  the  strength  of  the  optical  lattice  used  to  confine  the  gas,  in  units  of  the 
recoil  energy  ER  of  a  6Li  atom  in  the  lattice. 


2.  Revealing  the  Superfluid  Lambda  Transition  in  the  Universal  Thermodynamics 
of  a  Unitary  Fermi  Gas 

Mark  J.  FI.  Ku,  Ariel  T.  Sommer,  Lawrence  W.  Cheuk,  Martin  W.  Zwierlein 
Science  335,  563  (2012) 

Flighlighted  in  a  Science  Perspective  by  Wilhelm  Zwerqer 

Fermi  gases,  collections  of  fermions  such  as  neutrons  and  electrons,  are  found 
throughout  nature,  from  solids  to  neutron  stars.  Interacting  Fermi  gases  can  form  a 
superfluid  or,  for  charged  fermions,  a  superconductor.  We  have  directly  observed  the 
superfluid  phase  transition  in  a  strongly  interacting  Fermi  gas  via  high-precision 
measurements  of  the  local  compressibility,  density  and  pressure.  Our  data  completely 
determine  the  universal  thermodynamics  of  these  gases  without  any  fit  or  external 
thermometer.  The  onset  of  superfluidity  is  observed  in  the  compressibility,  the  chemical 
potential,  the  entropy,  and  the  heat  capacity,  which  displays  a  characteristic  lambda-like 
feature  at  the  critical  temperature  of  16.7%  of  the  Fermi  temperature.  Scaled  to  the 
density  of  electrons  in  a  metal,  this  form  of  superfluidity  would  occur  far  above  room 
temperature.  Our  measurements  provide  a  benchmark  for  many-body  theories  on  strongly 


interacting  fermions,  relevant  for  problems  ranging  from  high-temperature 
superconductivity  to  the  equation  of  state  of  neutron  stars. 


Figure  2  Observation  of  the  Superfluid  Lambda  Transition  in  a  trongly  interacting  Fermi  Gas.  Shown 
is  the  specific  heat  of  the  gas,  directly  obtained  from  the  density  profiles  of  a  trapped  gas. 

3.  Feynman  diagrams  versus  Feynman  quantum  emulator 

K.  Van  Floucke,  F.  Werner,  E.  Kozik,  N.  Prokofev,  B.  Svistunov,  M.  Ku,  A.  Sommer,  L.  W. 

Cheuk,  A.  Schirotzek,  M.  W.  Zwierlein 

Feynman  diagrams  versus  Fermi-gas  Feynman  emulator 

Nature  Physics  8,  366  (2012) 

Precise  understanding  of  strongly  interacting  fermions,  from  electrons  in  modern 
materials  to  nuclear  matter,  presents  a  major  goal  in  modern  physics.  Plowever,  the 
theoretical  description  of  interacting  Fermi  systems  is  usually  plagued  by  the  intricate 
quantum  statistics  at  play.  Plere  we  present  a  cross-validation  between  a  new  theoretical 
approach,  Bold  Diagrammatic  Monte  Carlo  (BDMC),  and  precision  experiments  on  ultra¬ 
cold  atoms.  Specifically,  we  compute  and  measure  with  unprecedented  accuracy  the 
normal-state  equation  of  state  of  the  unitary  gas,  a  prototypical  example  of  a  strongly 
correlated  fermionic  system.  Excellent  agreement  demonstrates  that  a  series  of  Feynman 
diagrams  can  be  controllably  resummed  in  a  non-perturbative  regime  using  BDMC.  This 
opens  the  door  to  the  solution  of  some  of  the  most  challenging  problems  across  many 
areas  of  physics. 


a)  b) 


Figure  3  Equation  of  State  of  a  Unitary  Fermi  Gas.  a)  Density  and  b)  pressure  as  a  function  of  the 
ratio  of  the  chemical  potential  to  the  temperature.  Density  and  pressure  are  normalized  by  the 
respective  quantities  of  a  non-interacting  Fermi  gas  at  the  same  ratio  p/kBT. 


4.  Quantum  degenerate  Bose-Fermi  mixture  of  chemically  different  atomic  species 
with  widely  tunable  interactions 

Jee  Woo  Park,  Cheng-Hsun  Wu,  Ibon  Santiago,  Tobias  G.  Tiecke,  Peyman  Ahmadi, 

Martin  W.  Zwierlein 

Phys.  Rev.  A  85,  051602(R)  (2012) 

We  have  created  a  quantum  degenerate  Bose-Fermi  mixture  of  23Na  and  40K  with  widely 
tunable  interactions  via  broad  interspecies  Feshbach  resonances.  Twenty  Feshbach 
resonances  between  23Na  and  40K  were  identified.  The  large  and  negative  triplet 
background  scattering  length  between  23Na  and  40K  causes  a  sharp  enhancement  of  the 
fermion  density  in  the  presence  of  a  Bose  condensate.  As  explained  via  the  asymptotic 
bound-state  model  (ABM),  this  strong  background  scattering  leads  to  a  series  of  wide 
Feshbach  resonances  observed  at  low  magnetic  fields.  Our  work  opens  up  the  prospect 
to  create  chemically  stable,  fermionic  ground  state  molecules  of  23Na-40K  where  strong, 
long-range  dipolar  interactions  will  set  the  dominant  energy  scale. 


(d)  23Na  (e)  «K 


Figure  4  Strongly  interacting  degenerate  Bose-Fermi  mixture  of  23Na  and  40K.  The  fermion  cloud  (to 
the  right)  displays  an  untypical  “bimodality”  due  to  the  strong  interactions  with  the  sodium 

condensate. 

5.  Spin-Injection  Spectroscopy  of  a  Spin-Orbit  Coupled  Fermi  Gas 

Lawrence  W.  Cheuk,  Ariel  T.  Sommer,  Zoran  Fladzibabic,  Tarik  Yefsah,  Waseem  S.  Bakr, 
Martin  W.  Zwierlein 

Phys.  Rev.  Lett.,  in  print,  preprint  arXiv:  1205.3483  (2012) 

The  coupling  of  the  spin  of  electrons  to  their  motional  state  lies  at  the  heart  of  recently 
discovered  topological  phases  of  matter.  Here  we  create  and  detect  spin-orbit  coupling  in 
an  atomic  Fermi  gas,  a  highly  controllable  form  of  quantum  degenerate  matter.  We  reveal 
the  spin-orbit  gap  via  spin-injection  spectroscopy,  which  characterizes  the  energy- 
momentum  dispersion  and  spin  composition  of  the  quantum  states.  For  energies  within 
the  spin-orbit  gap,  the  system  acts  as  a  spin  diode.  To  fully  inhibit  transport,  we  open  an 
additional  spin  gap,  thereby  creating  a  spin-orbit  coupled  lattice  whose  spinful  band 
structure  we  probe.  In  the  presence  of  s-wave  interactions,  such  systems  should  display 
induced  p-wave  pairing,  topological  superfluidity,  and  Majorana  edge  states. 


q/Q  q/Q  q/Q 

Figure  5  Spin-Injection  Spectroscopy  of  a  Spin-Orbit  Coupled  Fermi  Gas.  (A)  A  radiofrequency 
pulse  transfers  atoms  from  the  reservoir  states  |t>R  and  |I>R  (shown  in  black)  into  the  spin-orbit 
coupled  system  (shown  in  red  and  blue).  Transfer  occurs  when  the  RF  photon  energy  equals  the 
energy  difference  between  the  reservoir  state  and  the  spin-orbit  coupled  state  at  quasi-momentum 


q.  (B,C,D  and  E)  Transfer  as  a  function  of  RF  frequency  and  detuning  hAv  and  quasi-momentum  q. 
(B  and  C)  Spin-resolved  |T>  and  |1>  spectra,  respectively,  when  transferring  out  of  |t>R.  (D  and  E) 
Spin-resolved  |t>  and  |4>  spectra,  respectively,  when  transferring  out  of  |I>R.  (F,  G  and  FI)  The 
reconstructed  spinful  dispersions  of  a  spin-orbit  coupled  Fermi  gas  for  various  strengths  of  Raman 

coupling. 

6.  Ultracold  Fermionic  Feshbach  Molecules  of  23Na40K 

Cheng-Hsun  Wu,  Jee  Woo  Park,  Peyman  Ahmadi,  Sebastian  Will,  Martin  W.  Zwierlein. 
Phys.  Rev.  Lett.,  in  print,  preprint  arXiv:1 206.5023  (2012). 


We  report  on  the  formation  of  ultracold  fermionic  Feshbach  molecules  of  23Na40K,  the  first 
fermionic  molecule  that  is  chemically  stable  in  its  ground  state.  The  lifetime  of  the  nearly 
degenerate  molecular  gas  exceeds  100  ms  in  the  vicinity  of  the  Feshbach  resonance.  The 
measured  dependence  of  the  molecular  binding  energy  on  the  magnetic  field 
demonstrates  the  open-channel  character  of  the  molecules  over  a  wide  field  range  and 
implies  significant  singlet  admixture.  This  will  enable  efficient  transfer  into  the  singlet 
vibrational  ground  state,  resulting  in  a  stable  molecular  Fermi  gas  with  strong  dipolar 
interactions. 
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Figure  6  Radiofrequency  association  of  ultracold  Feshbach  molecules  of  NaK.  Starting  with  a 
mixture  of  sodium  atoms  in  hyperfine  state  |1,1>Na  and  potassium  atoms  in  hyperfine  state  |9/2,- 
3/2>k,  rf  spectroscopy  near  the  |9/2,-3/2>K  to  |9/2,-5/2>K  hyperfine  transition  reveals  free  40K  atoms 
repulsively  interacting  with  the  23Na  bath  (near  zero  rf  offset),  as  well  as  associated  molecules  (near 
85  kFIz  rf  offset).  A  fit  to  the  molecular  association  spectrum  yields  a  binding  energy  of  Eb  =  h  x 
84(6)  kHz.  The  magnetic  field  corresponding  to  the  atomic  transition  at  34.975  MHz  was  1 29.4  G. 
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Spin-Injection  Spectroscopy  of  a  Spin-Orbit  Coupled  Fermi  Gas 

Lawrence  W.  Cheuk,1  Ariel  T.  Sommer,1  Zoran  Hadzibabic,1, 2 
Tarik  Yefsah,1  Waseem  S.  Bakr,1  and  Martin  W.  Zwierlein1 
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2  Cavendish  Laboratory,  University  of  Cambridge,  J.  J.  Thomson  Avenue,  Cambridge  CB3  OHE,  United  Kingdom 

The  coupling  of  the  spin  of  electrons  to  their  motional  state  lies  at  the  heart  of  recently  discovered 
topological  phases  of  matter  [1-3].  Here  we  create  and  detect  spin-orbit  coupling  in  an  atomic 
Fermi  gas,  a  highly  controllable  form  of  quantum  degenerate  matter  [4,  5] .  We  reveal  the  spin-orbit 
gap  [6]  via  spin-injection  spectroscopy,  which  characterizes  the  energy-momentum  dispersion  and 
spin  composition  of  the  quantum  states.  For  energies  within  the  spin-orbit  gap,  the  system  acts  as  a 
spin  diode.  To  fully  inhibit  transport,  we  open  an  additional  spin  gap,  thereby  creating  a  spin-orbit 
coupled  lattice  [7]  whose  spinful  band  structure  we  probe.  In  the  presence  of  s-wave  interactions, 
such  systems  should  display  induced  p-wave  pairing  [8],  topological  superfluidity  [9],  and  Majorana 
edge  states  [10]. 


Spin-orbit  coupling  is  responsible  for  a  variety  of  phe¬ 
nomena,  from  the  fine  structure  of  atomic  spectra  to  the 
spin  Hall  effect,  topological  edge  states,  and  the  predicted 
phenomenon  of  topological  superconductivity  [3,  11].  In 
electronic  systems,  spin-orbit  coupling  arises  from  the  rel¬ 
ativistic  transformation  of  electric  fields  into  magnetic 
fields  in  a  moving  reference  frame.  In  the  reference 
frame  of  an  electron  moving  with  wavevector  k  in  an 
electric  field,  the  motional  magnetic  field  couples  to  the 
electron  spin  through  the  magnetic  dipole  interaction. 
This  spin-orbit  coupling  phenomenon  is  responsible  for 
lifting  the  degeneracy  of  spin  states  in  the  excited  or¬ 
bitals  of  atoms  and  solid-state  materials  such  as  zinc- 
blende  structures  [12].  In  a  two-dimensional  semicon¬ 
ductor  heterostructure,  the  electric  field  can  arise  from 
structure  or  bulk  inversion  asymmetry  [13],  leading  to 
magnetic  fields  of  the  form  =  a(—ky,kx,  0)  or 

B(£>)  =  (3(kv,kx,  0).  The  resulting  spin-orbit  coupling 
terms  in  the  Hamiltonian  are  known  as  the  Rashba  [14] 
and  Dresselhaus  [12]  contributions,  respectively.  In¬ 
cluding  a  possible  momentum-independent  Zeeman  held 
=  (0,  ByZ\  B^),  the  Hamiltonian  of  the  electron 
takes  the  form: 

U=— A  -  •  (B(£>) +B(fl) +B(Z)),  (1) 

2  m  li 

where  g  is  the  electron  g- factor,  is  the  Bohr  magneton 
and  S  is  the  electron  spin. 

The  energy-momentum  dispersion  and  the  associated 

spin  texture  of  the  Hamiltonian  in  Eq.  (1)  are  shown  in 
( Z) 

Figure  1A  for  By  —  0  and  a  =  ff.  In  the  absence  of  a 
perpendicular  Zeeman  held  B~  ,  the  spectrum  consists 
of  free  particle  parabolas  for  the  two  spin  states  that  are 
shifted  relative  to  each  other  in  fc-space  owing  to  the  spin- 
orbit  interaction.  For  a  hnite  held  BiZ\  a  gap  opens  in 
the  spectrum.  This  gap,  known  as  the  spin-orbit  gap, 
has  been  recently  observed  in  one-dimensional  quantum 
wires  [6,  15].  The  two  energy  bands  are  spinful  in  the 
sense  that  the  spin  of  an  atom  is  locked  to  its  momentum. 


Similar  band  structures  have  been  used  to  explain  the 
anomalous  quantum  Hall  effect  and  predict  a  saturation 
of  the  Hall  conductivity  for  Fermi  energies  in  the  gap 
region  [16]. 

In  this  work,  we  engineer  the  Hamiltonian  in  Eq.  (1) 
with  equal  Rashba  and  Dresselhaus  strengths  in  an  opti¬ 
cally  trapped,  degenerate  gas  of  fermionic  lithium  atoms 
via  Raman  dressing  of  atomic  hyperhne  states  [17,  18]. 
Raman  helds  have  previously  been  used  to  generate  spin- 
orbit  coupling  and  gauge  helds  in  pioneering  work  on 
Bose- Einstein  condensates  [19-21],  and  recently  spin- 
orbit  coupling  in  Fermi  gases  [22],  Here,  we  directly 
measure  the  spin-orbit  band  structure  of  Eq.  (1),  as  well 
as  the  rich  band  structure  of  a  spin-orbit  coupled  lat¬ 
tice.  For  this,  we  introduce  spin-injection  spectroscopy, 
which  is  capable  of  completely  characterizing  the  quan¬ 
tum  states  of  spin-orbit  coupled  fermions,  including  the 
energy-momentum  dispersion  and  the  associated  spin- 
texture.  By  tracing  the  evolution  of  quantum  states  in 
the  Brillouin  zone,  this  method  is  able  to  directly  mea¬ 
sure  topological  invariants,  such  as  the  Chern  number  in 
a  two-dimensional  system  [3,  11,  23]. 

In  order  to  directly  reveal  the  single-particle  eigen¬ 
states  of  the  spin-orbit  coupled  system,  we  reduce  the  in¬ 
teractions  in  our  Fermi  gas  to  a  negligible  strength.  This 
is  convenient  for  studying  topological  insulators,  whose 
behavior  is  mostly  governed  by  single-particle  physics. 
On  the  other  hand,  a  single-component  spin-orbit  cou¬ 
pled  Fermi  gas  is  expected  to  develop  effective  p-wave 
interactions  mediated  by  s-wave  interactions  [8],  either 
in  the  presence  of  an  s-wave  Feslrbaclr  resonance,  or 
in  the  presence  of  flat  bands  as  realized  below.  This 
can  lead  to  BCS  pairing  in  a  p-wave  channel,  and  in  a 
two-dimensional  system  with  pure  Rashba  coupling,  to 
Px  +  Wy  pairing  and  chiral  superfluidity  [8,  9]. 

We  generate  spin-orbit  coupling  using  a  pair  of  laser 
beams  that  connect  two  atomic  hyperfme  levels,  la¬ 
beled  |t)  and  |t),  via  a  two-photon  Raman  transition 
(Fig.  1B,C).  The  Raman  process  imparts  momentum 
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FIG.  1.  Realization  of  spin-orbit  coupling  in  an  atomic  Fermi  gas.  (A)  Energy  bands  as  a  function  of  quasi-momentum  q  for 
Raman  coupling  strength  of  =  0.25 Er  and  H5  =  0.  The  spin  composition  of  the  states  is  indicated  by  the  color.  Dashed 
lines  show  energy  bands  for  FiQr  =  0 Er  and  hS  =  0 Er.  (B)  Geometry  of  the  Raman  beams:  A  pair  of  Raman  beams  at  ±19° 
relative  to  the  y  axis  couples  states  |4,  kx  =  q)  and  \\,kx  =  q  +  Q) .  A  bias  magnetic  field  B  in  the  z  direction  provides  the 
quantization  axis.  (C)  The  hyperfine  interaction  splits  |'|")  and  |4)  by  ftuo,  and  the  relevant  polarization  components  are  7r  and 
HS  is  the  two-photon  detuning.  (D)  Momentum-dependent  Rabi  oscillations  in  the  spin  texture  after  sudden  switch-on  of 
the  Raman  beams.  Here  HOLr  =  0.78(2 )Er  and  the  detuning  HS  =  —0.25(1  )Er.  (E)  A  7r-pulse  for  the  resonant  momentum- 
class  of  atoms  is  applied  at  different  two-photon  detunings  hS.  The  Raman  strength  is  hO.R  =  0.035(5)Er  in  order  to  retain 
momentum  selectivity.  (F  and  G):  Adiabatic  loading  and  unloading  of  atoms  into  the  upper  (lower)  band  at  coupling  strength 
of  HQr  =  0.53(5)Er.  The  Raman  beams  are  turned  on  with  5  =  =f8.5B_r,  and  the  detuning  is  swept  linearly  to  <S  =  0  and  back 
at  a  rate  of  |j|  =  0.27(5)Bfj.  This  loads  atoms  into  the  upper  (lower)  band,  as  indicated  by  the  diagrams  on  the  right.  The 
spin  texture  follows  the  instantaneous  value  of  S,  indicating  adiabaticity. 


fiQx  to  an  atom  while  changing  its  spin  from  |4)  to  |f), 
and  momentum  —HQx  while  changing  the  spin  from  |f) 
to  |4,).  Defining  a  quasimomentum  q  =  kx  —  ^  for  spin 
|4.)  and  q  =  kx  +  ®  for  spin  |f),  one  obtains  the  Hamilto¬ 
nian  of  the  form  given  in  Eq.  (1)  [19].  In  this  mapping, 
Bz  '  =  HQji/gfiB,  where  Qr  is  the  two-photon  Rabi  fre- 
quency,  By  '  =  hS/ggs,  where  8  is  the  two-photon  de- 
tuning,  and  a  =  /3  =  2mgtB  (see  Supplemental  Material). 

When  the  spin-orbit  gap  is  opened  suddenly,  an  atom 
prepared  in  the  state  |4 ,kx  =  q  —  Q/2)  oscillates  between 
|4,  kx  =  q  —  Q/2)  and  |f,  kx  =  q  +  Q/2 )  with  a  momen¬ 
tum  dependent  frequency  A (q)/h,  where  A(g)  is  the  en¬ 
ergy  difference  between  the  bands  at  quasimomentum 
q.  Such  Rabi  oscillations  correspond  to  Larmor  preces¬ 
sion  of  the  pseudo-spin  in  the  effective  magnetic  field 
B(s°)  =  b(d)  +  B(r)  +  B^z).  We  have  observed  these 
oscillations  by  starting  with  atoms  in  |4),  pulsing  on  the 
Raman  field  for  a  variable  duration  r,  and  imaging  the 
atoms  spin-selectively  after  time-of-flight  expansion  from 
the  trap.  Time-of-flight  maps  momentum  to  real  space, 
allowing  direct  momentum  resolution  of  the  spin  popula¬ 
tions.  As  a  function  of  pulse  duration,  we  observe  oscil¬ 
lations  of  the  pseudospin  polarization  with  momentum- 
dependent  frequencies  (Fig.  ID).  Our  Fermi  gas  occupies 


a  large  range  of  momentum  states  with  near-unity  occu¬ 
pation.  Therefore,  each  image  at  a  given  pulse  duration  r 
contains  information  for  a  large  range  of  momenta  q.  The 
observation  of  momentum-dependent  oscillations  demon¬ 
strates  the  presence  of  a  spin-orbit  gap,  and  shows  that 
the  atomic  system  is  coherent  over  many  cycles.  To  high¬ 
light  the  momentum  selectivity  of  this  process,  we  pre¬ 
pare  an  equal  mixture  of  atoms  in  states  |f)  and  |4)  and 
pulse  on  the  Raman  fields  for  a  time  t  =  7r/f Ir  for  dif¬ 
ferent  two-photon  detunings  <5.  This  inverts  the  spin  for 
atoms  with  momentum  q  where  A (q)  is  minimal.  Since 
the  minimum  of  A  (q)  varies  linearly  with  S  due  to  the 
Doppler  shift  c x  kxQ,  the  momentum  q  at  which  the  spin 
is  inverted  depends  linearly  on  <5  (Fig.  IE). 

Instead  of  pulsing  on  the  Raman  field  and  projecting 
the  initial  state  into  a  superposition  of  states  in  the  two 
bands,  one  can  introduce  the  spin-orbit  gap  adiabatically 
with  respect  to  band  populations.  This  is  achieved  by 
starting  with  a  spin-polarized  Fermi  gas  and  sweeping 
the  two-photon  detuning  S  from  an  initial  value  Si  to  a 
final  detuning  Sf.  The  magnitude  of  the  initial  detun¬ 
ing  5,  is  much  larger  than  the  two-photon  recoil  energy 
Er  =  h2Q2 /2m,  so  that  the  effective  Zeeman  field  is  al¬ 
most  entirely  parallel  with  the  spins.  Depending  on  the 
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FIG.  2.  Spin-injection  spectroscopy  of  a  spin-orbit  coupled  Fermi  gas.  (A)  A  radiofrequency  (RF)  pulse  transfers  atoms  from 
the  reservoir  states  (shown  in  black)  1)*)^  and  |4-)H  into  the  spin-orbit  coupled  system  (shown  in  red  and  blue).  Transfer  occurs 
when  the  RF  photon  energy  equals  the  energy  difference  between  the  reservoir  state  and  the  spin-orbit  coupled  state  at  quasi¬ 
momentum  q.  (B,C,D  and  E)  Transfer  as  a  function  of  RF  frequency  detuning  hAv  and  quasi  momentum  q  at  Raman  coupling 
strength  of  =  0.43(5)15#  and  hS  =  0.00(3)i5#.  Note  that  starting  with  reservoir  |J.)H  (|t)fl)i  transfer  to  state  |1")  (If))  is 
entirely  due  to  spin-orbit  coupling.  Hence  the  signal  is  generally  much  weaker  than  that  for  state  |4-)  (|t))  except  right  in  the 
gap,  where  their  ratio  approaches  50%/50%.  (B  and  C)  Spin-resolved  |4-)  and  |i")  spectra,  respectively,  when  transferring  out  of 
(D  and  E)  Spin-resolved  |f)  and  |f)  spectra,  respectively,  when  transferring  out  of  If)^.  (F,G  and  H)  The  reconstructed 
spinful  dispersions  for  hS  =  0.00(3)15#  and  =  0 Er,  fifl#  =  0.43(5)15#  and  Ml#  =  0.9(l)i5#,  respectively. 


direction  of  the  sweep,  this  loads  atoms  into  either  the 
upper  or  the  lower  dressed  band.  We  interrupt  the  sweep 
at  various  times,  and  image  the  spin- momentum  distri¬ 
bution.  This  reveals  that  the  spin  texture  follows  the 
effective  Zeeman  field.  The  process  is  reversible,  as  we 
verify  by  sweeping  the  detuning  back  to  <5*  and  restoring 
full  spin-polarization.  (Fig.  IF  and  G). 

Having  demonstrated  the  ability  to  engineer  spin-orbit 
coupling  in  a  Fermi  gas,  we  introduce  a  general  approach 
to  measure  the  complete  eigenstates  and  energies  of 
fermions  at  each  quasi-momentum  q  and  thus  resolve  the 
band  structure  and  associated  spin  texture  of  spin-orbit 
coupled  atomic  systems.  Our  approach  yields  equivalent 
information  to  spin  and  angle-resolved  photoemission 
spectroscopy  (spin-ARPES),  a  powerful  technique  re¬ 
cently  developed  in  condensed  matter  physics  [24] .  Spin- 
ARPES  is  particularly  useful  for  studying  magnetic  and 
quantum  spin  Hall  materials;  it  has  been  used,  for  exam¬ 
ple,  to  directly  measure  topological  quantum  numbers  in 
the  Bq.ajSbj,  series,  revealing  the  presence  of  topological 
order  and  chiral  properties  [25]. 

Our  spectroscopic  technique  uses  radiofrequency  (RF) 


spin-injection  of  atoms  from  a  free  Fermi  gas  into  an 
empty  spin-orbit  coupled  system  using  photons  of  a 
known  energy  (Fig.  2A).  After  injection,  the  momentum 
and  spin  of  the  injected  atoms  are  analyzed  using  time  of 
flight  [26]  combined  with  spin-resolved  detection.  Atoms 
are  initially  loaded  into  one  of  two  free  “reservoir”  atomic 
states  |f)#  and  |f)#  which  can  be  coupled  to  the  states  |f) 
and  |t),  respectively,  via  the  RF  spin- injection  field,  with¬ 
out  changing  the  quasimomentum.  The  injection  occurs 
when  the  frequency  of  the  RF  pulse  matches  the  energy 
difference  between  the  spin-orbit  coupled  bands  and  the 
initial  reservoir  state  (see  Fig.  2A).  Spin-injection  from 
|t)#  (It)#)  populates  mostly  the  region  of  the  spin-orbit 
coupled  bands  with  a  strong  admixture  of  |f)  (It))  states. 
Thus,  the  use  of  two  reservoir  states  allows  us  to  measure 
both  the  |t)-rich  and  the  |t)-rich  parts  of  the  spin-orbit 
coupled  bands.  Following  the  injection  process,  the  Ra¬ 
man  beams  are  switched  off,  and  the  atoms  are  simulta¬ 
neously  released  from  the  trap.  After  a  sufficiently  long 
time  of  free  expansion,  the  density  distribution  gives  ac¬ 
cess  to  the  momentum  distribution,  which  we  measure 
using  state-selective  absorption  imaging.  By  counting 
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the  number  of  atoms  of  a  given  spin  and  momentum  as 
a  function  of  injection  energy,  we  determine  the  disper¬ 
sion  of  the  spin-orbit  coupled  bands  along  with  their  spin 
texture. 

The  topological  characteristics  of  the  bands,  which  are 
encoded  in  the  eigenstates,  can  be  extracted  from  the 
spin  and  momentum  composition.  For  our  spin-orbit  sys¬ 
tem  with  (5  =  0,  the  spin  of  the  eigenstates  is  confined 
to  the  y-z  plane  on  the  Bloch  sphere  because  the  effec¬ 
tive  magnetic  field  has  no  x  component.  More  general 
couplings  may  not  restrict  the  spin  to  a  great  circle  on 
the  Bloch  sphere,  in  which  case  at  least  two  spin  compo¬ 
nents  must  be  measured  for  a  complete  characterization 
of  the  bands.  This  can  be  achieved  by  rotating  the  dif¬ 
ferent  spin  components  onto  the  measurement  basis  with 
an  RF  pulse. 

Applying  spin-injection  spectroscopy,  we  have  mea¬ 
sured  the  band  structure  of  the  equal-part  Rashba- 
Dresselhaus  Hamiltonian  at  <5  =  0  for  several  12#.  Fig¬ 
ure  2B,  C,  D  and  E  show  spin-  and  momentum-  resolved 
spin-injection  spectra  obtained  with  atoms  starting  in  the 
|t)#  reservoir  (top  row)  and  starting  in  the  ||)H  reservoir 
(bottom  row),  for  the  case  M2#  =  0.43(5)E#  and  (5  =  0. 
The  (q,  f)  -H-  (—<7,  |)  symmetry  of  the  system  can  be  seen 
in  the  spectra  in  Fig.  2.  The  energy  at  each  quasimo¬ 
mentum  is  found  by  adding  the  energy  injected  into  the 
system  by  the  RF  pulse  to  the  initial  kinetic  energy  of  the 
free  particle  in  the  reservoir.  Figure  2F,  G  and  H  show 
the  dispersion  and  spin  texture  of  the  bands  obtained 
from  the  data.  As  f 1#  is  increased,  we  observe  the  open¬ 
ing  of  a  spin-orbit  gap  at  q  =  0.  The  spin  composition  of 
the  bands  evolves  from  purely  |f)  or  ||)  away  from  the 
spin-orbit  gap  to  a  mixture  of  the  two  spin  states  in  the 
vicinity  of  the  spin-orbit  gap,  where  the  spin  states  are 
resonantly  coupled. 

The  dispersion  investigated  above  is  the  simplest  pos¬ 
sible  for  a  spin-orbit  coupled  system  and  arises  naturally 
in  some  condensed  matter  systems.  A  Fermi  gas  with 
this  dispersion  has  an  interesting  spinful  semi-metallic 
behavior  when  the  Fermi  energy  lies  within  the  spin- 
orbit  gap.  When  the  Fermi  energy  is  outside  the  spin- 
orbit  gap,  there  is  a  four-fold  degeneracy  of  states  at  the 
Fermi  surface.  Inside  the  gap,  however,  the  degeneracy 
is  halved.  Furthermore,  propagation  of  spin  up  parti¬ 
cles  at  the  Fermi  energy  can  only  occur  in  the  positive 
q  direction,  while  spin  down  fermions  can  only  propa¬ 
gate  in  the  opposite  way.  Particles  are  thus  protected 
from  back-scattering  in  the  absence  of  magnetic  impuri¬ 
ties  that  would  rotate  their  spin.  Such  a  spinful  semi¬ 
metal  can  be  used  to  build  spin-current  diodes,  since  the 
material  permits  flow  of  polarized  spin-currents  in  one 
direction  only. 

An  even  richer  band  structure  involving  multiple  spin¬ 
ful  bands  separated  by  fully  insulating  gaps  can  arise 
in  the  presence  of  a  periodic  lattice  potential.  This  has 
been  realized  for  Bose-Einstein  condensates  by  adding 


RF  coupling  between  the  Raman-coupled  states  |f)  and 
||)  [7].  Using  a  similar  method,  we  create  a  spinful  lat¬ 
tice  for  ultracold  fermions,  and  use  spin-injection  spec¬ 
troscopy  to  probe  the  resulting  spinful  band  structure. 
The  combined  Raman/RF  coupling  scheme  is  shown  in 
Fig.  3A.  The  Raman  field  couples  the  states  ||,  kx  =  q ) 
and  |t ,kx  =  q  +  Q)  with  strength  f 1#,  whereas  the  RF 
field  couples  the  states  ||,  kx  =  q)  and  |f,  kx  =  q)  with 
strength  12##.  As  a  result,  the  set  of  coupled  states  for 
a  given  quasimomentum  q,  shown  in  the  repeated  Bril- 
louin  scheme  in  Fig.  3B,  is  |cr,  kx  =  q  +  nQ)  for  integer 
n  and  a  =f,|.  The  lowest  four  bands  are  degenerate  at 
the  band  center  q  =  0  when  12#  =  12##  =  0.  The  Raman 
field  splits  the  degeneracy  between  the  first  and  fourth 
band,  leaving  the  other  two  degenerate.  The  remaining 
degeneracy,  which  is  a  Dirac  point,  is  removed  with  the 
addition  of  the  RF  field.  Thus,  when  the  system  is  filled 
up  to  the  top  of  the  second  band,  it  is  an  insulator.  Fur¬ 
thermore,  when  12##  is  large  enough,  a  band  gap  also 
opens  between  the  first  and  second  bands. 

Fig.  3D  and  E  show  the  spin-injection  spectra,  mea¬ 
sured  with  fermions  initially  in  reservoir  state  ||)#,  which 
is  sufficient  to  reconstruct  the  full  band  structure  given 
the  (g,  t)  -H-  (—<7,1)  symmetry  of  the  Hamiltonian.  The 
transitions  between  the  reservoir  and  the  spin-orbit  cou¬ 
pled  bands  for  M2#  =  0.40(5)E#  and  Ml =  0.28(2 )ER 
are  shown  in  Fig.  3C.  The  experimental  spectra  (Fig.  3D 
and  E)  for  the  same  parameters  are  compared  to  the 
corresponding  theoretically  calculated  spectra,  shown  in 
Fig.  3F  and  G.  We  focus  on  the  features  of  the  ||)  channel 
of  the  spectrum,  which  is  stronger  because  of  the  better 
spin-composition  overlap  with  the  reservoir  state.  The 
spectrum  exhibits  four  prominent  features  separated  by 
three  energy  gaps,  labeled  Ai,  A2  and  A3  in  Fig.  3F  and 
3G.  The  gaps  giving  rise  to  these  features  are  shown  on 
the  band  structure  in  Fig.  3C.  The  gap  Ai  is  opened  by 
the  spin-orbit  coupling,  while  A2  is  opened  by  a  direct 
RF  coupling  and  A3  is  opened  by  a  second  order  process 
that  involves  both  the  RF  and  Raman  fields,  explaining 
its  smallness.  We  have  explored  the  Raman/RF  system 
for  a  range  of  coupling  strengths  as  shown  in  the  spectra 
in  Fig.  4B  and  4C.  The  corresponding  band  structures  are 
shown  in  Fig  4A.  With  a  careful  choice  of  the  Raman/RF 
coupling  strengths,  spinful  flat  bands  are  realized,  where 
interactions  should  play  a  dominant  role  [27]. 

To  illustrate  how  the  energy  bands  along  with  the 
corresponding  eigenstates  can  be  extracted,  we  recon¬ 
struct  the  energy  bands  along  with  the  spin  texture  for 
M2#  =  0.93(7)E#  and  M2##  =  0.28(2)1/#,  as  shown  in 
Fig.  4D.  The  energies  of  the  bands  are  obtained  from 
the  resonant  frequencies  in  the  spin-injection  spectra, 
while  the  spin  composition  is  extracted  from  the  relative 
weights  of  the  signal  in  the  two  spin  channels  (see  Supple¬ 
mental  Material).  In  general,  the  eigenvector  \4>^n\q)) 
for  the  nth  energy  band  at  a  given  quasi-momentum  q , 
can  be  expanded  in  terms  of  free  space  eigenstates  as 
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FIG.  3.  Creating  and  probing  a  spin-orbit  coupled  lattice.  (A)  The  addition  of  a  radiofrequency  field  allows  momentum  transfer 
of  any  multiple  of  Q.  The  combined  Raman-RF  system  produces  a  spinful  lattice  band  structure.  (B)  The  band  structure  of 
the  Raman-RF  system  in  the  repeated  zone  scheme.  The  topmost  band  structure  corresponds  to  hQ.RF  =  0  and  HQr  =  0.25 Er, 
which  has  a  band  crossing  at  quasi-momentum  q  =  0.  The  middle  band  structure  corresponds  to  a  larger  Raman  coupling  of 
HQr  =  0.5 Er  with  HQrf  =  0.  In  the  bottom-most  band  structure,  HQr  =  0.5Er  while  HQrf  is  increased  to  0.25Er.  (C)  Spin 
injection  from  free  particle  bands  to  spinful  lattice  bands,  starting  from  14,}^.  Transitions  near  zero  RF  detuning  (hAv  ~  0) 
that  give  rise  to  dominant  spectral  features  are  identified.  (D  and  E)  Experimental  spectrum  of  the  Raman-RF  system  with 
hflR  =  0.40(5).Er  and  FiQrf  =  0.28(2 )Er  in  the  spin  |1)  and  spin  |'j")  channels,  measured  after  injection  from  reservoir 
The  dominant  features  span  many  Brillouin  zones,  corresponding  to  projection  of  lattice  states  onto  free  particle  states  after 
time-of- flight.  (F  and  G)  The  theoretical  spectra  corresponding  to  D  and  E,  respectively.  The  features  corresponding  to  the 
gaps  and  transitions  identified  in  C  are  labeled. 


I ipn(q)}  =  Em, <7  c^\kx  =  <7  +  m Q )  |cr,  kx  =  q  +  mQ) .  In 
spin-injection  spectroscopy,  the  projection  of  the  lattice 
wavefunctions  onto  free  particle  states  allows  us  to  not 
only  extract  the  average  spin,  but  also  the  magnitude  of 
the  coefficients  c^l\kx).  From  the  projection  coefficients 
Co^ikx),  one  can  define  the  spin  S(kx)  (see  Supplemental 
materials).  In  Fig.  4E,  F  and  G,  we  show  the  extracted 
value  of  Sy{kx )  and  |5'2:(fc2,)|  for  the  bottommost  band 
when  =  0.93(7)7?^  and  MIrf  =  0.28(2 )Er.  For 
more  general  spin-orbit  Hamiltonians  involving  ax,  one 
can  extract  the  phase  between  all  three  components  of 
S(kx)  with  additional  RF  pulses,  and  fully  characterize 
the  eigenstate  for  the  corresponding  quasimomentum  q. 
The  topology  of  the  band,  encoded  in  the  evolution  of  its 
eigenstates  across  the  Brillouin  zone,  can  thus  be  mea¬ 
sured. 

In  summary,  we  have  created  and  directly  probed  a 
spin-orbit  gap  in  a  Fermi  gas  of  ultracold  atoms  and  re¬ 
alized  a  fully  gapped  band  structure  allowing  for  spinful 
flat  bands.  We  introduced  spin-injection  spectroscopy 
to  characterize  the  spin-textured  energy-momentum  dis¬ 
persion.  Such  measurements  would  reveal  the  non-trivial 
topology  of  the  bands  in  systems  with  more  general  spin- 
orbit  couplings  [28],  opening  a  path  to  probing  topo¬ 
logical  insulators  with  ultracold  atoms.  Recently  de¬ 


veloped  high  numerical  aperture  imaging  techniques  can 
be  used  for  microscopic  patterning  of  lower  dimensional 
Fermi  gases  into  heterostructures  with  regions  charac¬ 
terized  by  different  topological  numbers  separated  by 
sharp  interfaces  [29,  30].  In  such  systems,  spatially 
resolved  spin-injection  spectroscopy  can  directly  reveal 
topologically  protected  edge  states  such  as  Majorana 
fermions,  which  have  been  proposed  for  topological  quan¬ 
tum  computation[10,  31,  32]. 
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Army  Research  Office  with  funding  from  the  DARPA 
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MURI,  ONR  YIP,  DARPA  YFA,  an  AFOSR  PECASE, 
and  the  David  and  Lucile  Packard  Foundation.  Z.  H. 
acknowledges  funding  from  EPSRC  under  Grant  No. 
EP/I010580/1. 
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FIG.  4.  Evolution  of  spin-textured  energy  bands  of  a  spin-orbit  coupled  lattice.  (A)  Theoretical  band  structures  for  various 
combinations  of  Qr  and  Qrf-  The  first  band  becomes  flat  while  remaining  spinful  for  hflR  =  0.93 Er  and  H^Irf  =  Q.IIEr  and 
0.28-E.r.  (B  and  C)  The  corresponding  experimental  Raman-RF  spin-injection  spectra  for  injection  from  for  channels  |4-) 
and  |t),  respectively.  The  color  map  used  is  the  same  as  Fig.  2B  and  2E  after  rescaling  to  the  maximum  intensity  (which  for 
4C  is  20%  of  4E),  except  for  the  top  left  panel  in  (C),  which  is  scaled  to  the  maximum  intensity  of  the  corresponding  panel  in 
(B).  Possible  interaction  effects  between  |t)  with  |t)fl  (see  Fig.  S2)  makes  only  the  dominant  features  resolvable  in  |t),  while 
finer  features  are  visible  in  |4-}.  (D)  Experimentally  reconstructed  band  structure  for  hQ.R  =  0.93(7 )Er  and  hQ.RF  =  0.28(2). 
The  spin  texture  is  indicated  by  the  color  of  the  points.  (E,F  and  G)  Experimentally  measured  spin  components  Sy  and  |57| 
as  a  function  of  momentum  kx  for  the  lattice  wavefunctions  corresponding  to  the  bottommost  band  in  D. 
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SUPPLEMENTAL  MATERIALS 
System  preparation  and  Raman  setup 

Fermionic  6Li  in  the  hyperfine  state 
|F  =  3/2,  rriF  =  3/2)  is  sympathetically  cooled  to 
degeneracy  by  23Na  atoms  in  a  magnetic  trap.  The 
atoms  are  transferred  to  a  nearly  spherical  crossed  opti¬ 
cal  dipole  trap  with  mean  trapping  frequency  ~  150  Hz. 
Depending  on  the  measurement,  the  atoms  are  then 
transferred  into  one  of  the  lowest  four  hyperfine  states 
in  the  ground  state  manifold  via  radiofrequency  (RF) 
sweeps.  The  four  lowest  hyperfine  states  are  states 
|j,),  It))  lt)_R  in  the  text.  Subsequently,  the  magnetic 
field  is  ramped  to  B  =  BqZ,  with  Bq  =  11.6  G.  The 
geometry  of  the  Raman  beams  is  shown  in  Fig.  IB.  The 
two  beams,  detuned  by  3.96  GHz  to  the  blue  of  the  Di 
line,  generate  a  moving  lattice  with  lattice  wavevector 
Q  =  2tt  x  (1.0/nn)-1  and  corresponding  recoil  energy 
of  Er  =  — =  h  x  32(1)  kHz.  The  calibration  of 
the  recoil  energy  Er  is  performed  using  the  data  in 
Fig.  IE  and  relies  only  on  the  relation  of  Q  to  the 
traveled  distance  for  a  given  time-of- flight.  To  couple 
states  |t)  and  |t),  the  frequency  difference  between  the 
two  beams  is  set  near  the  hyperfine  splitting  of  ojq  = 
2i rx  207.7  MHz.  For  the  one-photon  detuning  that  we 
use,  0^/rsc  ss  240,  where  Tsc  is  the  scattering  rate. 
Note  that  one-body  losses  from  single-photon  scattering 
events  do  not  perturb  the  measured  spin-injection 
spectra,  as  they  affect  all  energy  states  equally  and  only 
reduce  the  number  of  atoms  in  each  momentum  state. 


Experimental  procedure 

In  the  presence  of  the  Raman  beams,  a  differential 
Stark  shift  between  |f)  and  |/)  can  alter  the  resonant 
two-photon  frequency.  The  frequency  corresponding  to 
6  =  0  is  calibrated  using  RF  spectroscopy  on  the  |f)  to 
|j,)  transition  in  the  presence  of  the  Raman  beams  with 
a  large  two-photon  detuning  8  «  27rxl  MHz.  The  res¬ 
onant  frequencies  for  Ij.)^  — »  |j,)  and  |t)^  — >  |t)  are 
calibrated  similarly.  The  differential  Stark  shift  for  the 
largest  Raman  coupling  strength  is  measured  to  be  4  kHz 
for  |t)  — >  |t))  and  <  1  kHz  for  the  |t)#  — >  |t)  and 
lt)*->  |t)  transitions.  For  spin-injection  spectroscopy, 
the  injection  process  uses  a  RF  field  that  couples  the 
states  |t)fl  (It)#)  and  |t)  (|j.)fi).  All  spectra  are  taken 
with  a  RF  injection  pulse  duration  of  0.5  ms,  and  an  in¬ 
jection  field  strength  corresponding  to  a  maximum  trans¬ 
fer  fraction  <  0.30.  The  experimental  frequency  resolu¬ 
tion  for  spin-injection  is  3  kHz  «  O.IEr,  while  the  mo¬ 
mentum  resolution  is  estimated  to  be  0.05Q,  limited  by 
expansion  time  and  imaging  resolution.  For  all  measure¬ 
ments,  state-selective  absorption  images  are  taken  after 


time-of-flight  of  4  ms.  As  the  atoms  are  released,  the 
bias  magnetic  field  is  ramped  to  300  G,  where  the  reso¬ 
nant  imaging  frequencies  for  different  hyperfine  states  are 
well-resolved,  allowing  spin-selective  absorption  imaging. 
To  obtain  the  spectra  shown  in  Fig.  1  and  2,  the  time-of- 
flight  images  for  each  spin  state  are  first  integrated  along 
y ,  orthogonal  to  the  spin-orbit  direction.  For  a  given 
quasi-momentum  q,  the  integrated  one-dimensional  den¬ 
sity  profiles  from  the  two  spin  channels  are  then  combined 
to  produce  the  final  spectrum.  As  an  example,  we  show 
in  Fig.  SI  the  time-of-flight  images  and  the  corresponding 
integrated  density  profiles,  at  a  specific  detuning  S/Er 
for  the  spectrum  in  Fig.  IE. 


Hamiltonian  for  Raman-coupled  system 


The  spin  of  an  atom  is  coupled  to  its  momentum  using 
a  pair  of  laser  beams  near  a  Raman  transition.  Using  the 
rotating  wave  approximation,  the  Raman  beams  generate 
a  spinor  potential 

HCIr  8 

V(r)  =  — —  (crx  cos  Qx  -  ay  smQx)  +  -<tz,  (2) 

where  Qr  is  the  two-photon  Rabi  frequency.  After  a  local 
pseudo-spin  rotation  about  the  z  axis  with  angle  Qx,  the 
Hamiltonian  becomes  [19] 


_  h2 k2  ,  h2Q  ,  hfln  s  ,  er 
rL  —  — - h  "x — crzq  H  —  &x  +  2°'z  4 — >  (3) 


2m 


2  m 


where  q  is  the  quasi-momentum  defined  in  the  text.  Fol¬ 
lowing  a  global  pseudo-spin  rotation  az  — >  ay,  <ry  — t  ax 
and  cjx  — >  az,  the  Hamiltonian  becomes 


n  = 


h2  k2 
2  m 


h2Q 
2  m 


<xyq  + 


~Y~ 


(4) 


which  up  to  a  constant  has  the  same  Rashba-Dresselhaus 
form  as  Eq.  (1),  with  a  =  /3  =  2^B ,  B(yZ)  =  hS/gyB 

and  Bz  =  Hflu/gfis-  In  this  convention,  the  bare  hy¬ 
perfine  states,  labeled  |f)  and  |j,)  are  eigenstates  of  ay. 


Reconstructing  the  spinful  dispersion  for  the 
Raman-coupled  system 

For  the  equal  Rashba-Dresselhaus  system,  the  two 
channels  in  a  spectrum  from  reservoir  er)  R  are  first  re¬ 
labeled  by  quasi-momentum  q.  The  ratio  of  the  trans¬ 
ferred  atoms  in  each  channel  at  a  given  q  directly  mea¬ 
sures  the  ^-dependent  spin  composition.  The  dispersion 
is  then  reconstructed  by  adding  the  free  particle  disper¬ 
sion  e°(q)  =  (q  =F  Q/ 2)2/2m  to  the  spectrum  correspond¬ 
ing  to  injection  from  |t),R  (|4-)ij)-  The  result  is  shown  in 
Fig.  2F,G  and  H,  where  the  color  denotes  the  spin  tex¬ 
ture  and  the  strength  of  the  color  is  weighed  by  the  total 
number  of  atoms  at  a  given  q. 
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FIG.  SI.  Example  of  converting  time-of-flight  images  to  a  spinful  spectrum.  Here,  we  show  the  conversion  process  for  the 
spectrum  in  Fig.  IE  for  a  specific  detuning  h5  =  —1.88 ER.  (A  and  B)  Time-of-flight  images  rescaled  in  terms  of  the  recoil 
momentum  Q  in  the  |f)  and  |4)  channels,  respectively.  (C  and  D)  The  time-of-flight  images  in  A  and  B  are  integrated  along  the 
ky  direction  to  produce  1-dimensional  densities  n(q/Q).  The  conversion  to  quasi-momentum  q  involves  adding  a  spin  channel 
dependent  momentum  offset.  (E)  Combining  C  and  D  for  every  quasi-momentum  q  produces  a  slice  in  the  final  spectrum.  In 
the  figure,  this  is  the  region  bounded  by  the  two  dashed  lines.  Repeating  the  same  procedure  for  other  detunings  yields  the 
full  spectrum. 


Hamiltonian  for  the  Raman/RF  system 


Spin-injection  spectrum  for  the  Raman/RF  system 


Raman  dressing  creates  a  spin-orbit  gap  in  momentum 
space.  Adding  radiofrequency  (RF)  coupling  with  zero 
momentum  transfer  creates  a  lattice  potential  with  true 
band  gaps  [7].  The  RF  drive  is  applied  at  the  same  fre¬ 
quency  as  the  Raman  frequency,  with  coupling  strength 
f Irf-  An  atom  can  now  receive  an  arbitrary  number  of 
units  of  the  recoil  momentum  HQ  by  interacting  alter¬ 
nately  with  the  Raman  field  and  the  RF  field  (see  Figure 
3(A)).  When  the  two  spin  states  are  coupled  at  different 
momenta  using  Raman  lasers,  and  additionally  at  the 
same  momentum  using  RF,  the  Hamiltonian  becomes 


h2k 2  HQr  ■  \  ,  mRF  ,  S 

ft  =  -^—  +  —^—{axcosQx-<TysmQx)-\ — c tx+-<tz . 


(5) 

Since  the  Hamiltonian  has  discrete  translational  symme¬ 
try  along  x,  its  eigenstates  can  be  expanded  in  plane 
waves  as 


|  'lftn(kx')}  — 


E  < 

■?>=Td 


.(*») 


(kx) 


o',  kx  =  kx  +jQ 


(6) 


where  kx  is  the  quasi-momentum  given  by  to  kx  restricted 
to  the  first  Brillouin  zone,  and  n  is  the  band  index. 


Starting  from  reservoir  er  =f,  4-,  the  injected  popula¬ 
tion  in  band  n  with  quasi-momentum  kx  at  a  given  RF 
frequency  w,  Pi  a{io,n,m,kx ),  is  given  by 

n,  l,  kx)  oc  n]a^2,na(kx  +  lQ)\c%l)(kx  +  jQ)\2 

3 

xC((hui  +  e°  (kx  +  jQ ))  —  (hojo  +  en(kxjj!}) 

where  f \i  G  is  the  RF  strength  coupling  the  reservoir  state 
to  |<t),  na{k)  is  the  trap-averaged  momentum  distribution 
for  reservoir  state  \<j) R,  tf(k)  =  h  t^ie  ^ree  Par“ 

tide  dispersion  of  the  /th  non-interacting  band,  en(k)  is 
the  dispersion  for  the  nth  band,  ujq  is  the  hyperfine  fre¬ 
quency  difference  between  |f)  and  ||),  and  C(x)  is  the  RF 
lineshape.  After  injection,  the  atoms  are  released  from 
the  trap.  After  sufficient  time-of-flight,  the  momentum 
distribution  is  given  by  the  real  space  atomic  density  pro¬ 
file,  which  for  the  spin  a'  channel  is 

PTOF,*',a(U’kx)  =  ^Pl,a(u,n3,kX)\c<£)  (kx)\2.  (8) 

n,l 

The  theoretical  spectra  in  Fig.  3F  and  G  are  obtained 
using  Eq.  (8)  and  coefficients  Ca(kx)  found  by  numeri¬ 
cally  diagonalizing  the  Hamiltonian  in  Eq.  (5).  The  ex¬ 
perimental  and  theoretical  spectra  for  Fig.  4  are  shown 
in  Fig.  S2. 
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FIG.  S2.  Experimental  and  theoretical  Spin-injection  spectra  of  Raman/RF  system  for  different  Raman/RF  strengths.  (A  and 
B)  The  experimental  Raman/RF  spin-injection  spectra  for  injection  from  IJ.)^  for  channels  |j_)  and  |j"),  respectively.  The  color 
map  used  is  the  same  as  Fig.  2B  and  2E  after  rescaling  to  the  maximum  intensity  (which  for  4C  is  20%  of  4E),  except  for 
the  top  left  panel  in  (B),  which  is  scaled  to  the  maximum  intensity  of  the  corresponding  panel  in  (A).  (C  and  D)  Theoretical 
spectra  corresponding  to  B  and  C.  We  do  not  take  into  account  finite  imaging  resolution,  which  affects  the  sharpness  along  the 
momentum  axis  kx/Q.  Scattering  due  to  residual  interactions  between  hyperfine  states  can  also  lead  to  blurring  of  momenta 
along  kx/Q ,  and  will  be  most  pronounced  for  the  spectra  in  B,  as  state  |j'}  scatters  more  strongly  with  14-)^  (scattering  length 
~  — 450ao)  than  |4_)  does  with  IJ.)^  (scattering  length  ~  — 2oo  with  typical  initial  1  /fez?  ~  3000ao  before  TOF). 


Reconstructing  the  spinful  band  structure  for  the 
Raman/RF  System 

We  first  describe  a  general  procedure  to  reconstruct  the 
band  structure  for  any  spinful  lattice  from  spin-injection 
spectra.  One  performs  spin-injection  spectroscopy  with 
reservoir  states  that  are  filled  up  to  a  Fermi  momentum 
smaller  than  half  the  recoil  momentum,  Ef  <  Q / 2.  One 
then  selects  a  prominent  feature  on  the  spectrum  and 


finds  the  resonant  transfer  frequencies  as  a  function  of  kx . 
over  one  Brillouin  zone.  A  mask  centered  on  the  resonant 
frequencies  and  repeated  over  all  kx  with  period  Q  is  then 
created.  With  the  mask  for  a  specific  feature  applied,  the 
spin  channel  a'  for  the  spectrum  starting  from  |  cr)  R  has 
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transfer  intensity  given  by 


I*'An,l,kx)  oc  Y,Mkx+lQ)\c^\kx+jQ)\2 3 4 5\c^(kx)\2. 
j 

(9) 

Here,  the  Fermi  momentum  kp  is  less  than  half  the  recoil 
momentum  Q / 2  and  therefore  only  the  l  =  0  contributes. 
Defining  Afa(n,  l,  kx)  =  A'An->  k  kx+jQ),  one  ob¬ 
tains 


c{:\kx)? 


k fj f  ,(T  {jl-,  l ,  kx) 

N<7  (jl,  k  kx) 


(10) 


Defining  S(kx)  =  \An\kx)Ac^n\kx),  where  c^n\kx)  = 
(c^\kx),  c^l\kx))T,  allows  to  measure  Sz(kx).  With  ad¬ 
ditional  RF  pulses,  one  can  measure  Sx(kx)  and  Sy(kx). 
After  the  bands  that  coupled  strongly  to  |cr,  k) ,  k  <  Q  are 
measured,  one  iterates  the  process  with  a  larger  Fermi  sea 
to  obtain  other  bands. 

In  the  text,  we  demonstrate  extraction  of  band  struc¬ 
ture  and  spin  texture  in  a  spin-orbit  coupled  lattice.  The 
initial  Fermi  sea  has  /cf  >  Q,  therefore  spectral  fea¬ 
tures  corresponding  to  transitions  out  of  different  non¬ 
interacting  reservoir  bands  can  all  appear  near  zero  RF 
detuning.  One  can  however  identify  spectral  features  cor¬ 
responding  to  different  transitions  and  apply  the  above 
procedure.  In  addition,  since  the  Raman/RF  system  is 
set  to  5  =  0,  one  can  invoke  the  additional  symmetry 
of  the  band  structure  about  q  =  0,  meaning  that  it  is 
sufficient  to  fit  a  certain  spectral  feature  over  half  of  a 
Brillouin  zone.  Due  to  a  larger  signal  in  the  |j,)  channel, 
features  in  this  channel  were  used  to  create  the  mask. 

In  order  to  be  consistent  with  the  earlier  conven¬ 
tion  for  the  Raslrba-Dresselhaus  system,  we  apply  a 
global  spin  rotation  such  that  Sy(kx)  corresponds  to 

\  ^|cjn'>(fca;)|2  —  |c[n\fca;)|2^ .  The  experimentally  mea¬ 
sured  Sy{kx)  for  the  lowest  band  for  MI#  =  0.93(7 )Ep 
and  hflpF  =  0.28(2 )Er  is  shown  in  Fig.  4F.  Since 
the  Hamiltonian  only  has  components  in  the  pseudo¬ 
spin  y  and  z  direction,  one  can  also  extract  |  •S's  ( )  |  = 

A/  I  —  Sy(kx)2,  as  shown  in  Fig.  4E. 
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will  slow  down  their  spin  rates  with  time  after  the 
RLDP.  However,  they  cannot  explain  the  apparent 
difference  in  spin  distributions  between  AXMSPs 
and  radio  MSPs,  because  radio  MSPs,  which 
have  weak  surface  magnetic  field  strengths,  could 
not  spin  down  by  the  required  amount  even  in  a 
Hubble  time.  The  true  age  of  a  pulsar  (23)  is  given 
by  t  =  P/{(n  -  1).P)[1  -  ( Po/P)n  *].  Assum¬ 
ing  an  evolution  with  a  braking  index  n  =  3 
and  B  =  l.Ox  10s  G,  the  time  scale  t  is  larger 
than  10  Gy,  using  P0  =  (P)axmsp  =  3.3  ms  and 
P(t)  =  (P) msp  =  5.5  ms.  To  make  things  worse, 
one  has  to  add  the  main-sequence  lifetime  of  the 
LMXB  donor  star,  which  is  typically  3  to  12  Gy, 
thereby  reaching  unrealistic  large  total  ages.  Al¬ 
though  the  statistics  of  AXMSPs  still  has  its  basis 
in  small  numbers  and  care  must  be  taken  for  both 
detection  biases  (such  as  eclipsing  effects  of  radio 
MSPs)  and  comparison  between  various  sub¬ 
populations  (5),  it  is  evident  from  both  observa¬ 
tions  and  theoretical  work  that  the  RLDP  effect 
presented  here  plays  an  important  role  for  the 
spin  distribution  of  MSPs. 

The  RLDP  effect  may  also  help  explain  a  few 
other  puzzles,  for  example,  why  characteristic  (or 
spin-down)  ages  of  radio  MSPs  often  largely  ex¬ 
ceed  cooling  age  determinations  of  their  white 
dwarf  companions  (24).  It  has  been  suggested 
that  standard  cooling  models  of  white  dwarfs  may 
not  be  correct  (25-27),  particularly  for  low-mass 
helium  white  dwarfs.  These  white  dwarfs  avoid 
hydrogen  shell  flashes  at  early  stages  and  retain 
thick  hydrogen  envelopes,  at  the  bottom  of  which 
residual  hydrogen  burning  can  continue  for  sev¬ 
eral  billion  years  after  their  fonnation,  keeping  the 
white  dwarfs  relatively  hot  (~104  K)  and  thereby 
appearing  much  younger  than  they  actually  are. 
However,  it  is  well  known  that  the  characteristic 
age  is  not  a  trustworthy  measure  of  true  age  (28), 
and  the  RLDP  effect  exacerbates  this  discrepancy 
even  further.  In  the  model  calculation  presented 
in  Fig.  1,  it  was  assumed  that  B  =  1.0  x  108G 
and  (p  =  1.0.  However,  P0  and  To  depend  strong¬ 
ly  on  bothf?  and  (p.  This  is  shown  in  Fig.  2,  where 
I  have  calculated  the  RLDP  effect  for  different 
choices  of  B  and  cp  by  using  the  same  stellar 
donor  model  [i.e.,  same  M(t)  profile]  as  before. 
The  use  of  other  LMXB  donor  star  masses,  met- 
allicities,  and  initial  orbital  periods  would  lead 
to  other M(t)  profiles  (16, 17)  and  hence  different 
evolutionary  tracks.  The  conclusion  is  that  recycled 
MSPs  can  basically  be  bom  with  any  characteristic 
age.  Thus,  we  are  left  with  the  cooling  age  of  the 
white  dwarf  companion  as  the  sole  reliable,  although 
still  not  accurate,  measure  as  an  age  indicator. 

A  final  puzzle  is  why  no  sub-millisecond  pul¬ 
sars  have  been  found  among  the  216  radio  MSPs 
detected  in  total  so  far.  Although  modem  obser¬ 
vational  techniques  are  sensitive  enough  to  pick 
up  sub-millisecond  radio  pulsations,  the  fastest 
spinning  known  radio  MSP,  J1748-2446ad  (29), 
has  a  spin  frequency  of  only  716  Hz,  correspond¬ 
ing  to  a  spin  period  of  1.4  ms.  This  spin  rate  is 
far  from  the  expected  minimum  equilibrium  spin 
period  (8)  and  the  physical  mass  shedding  limit 


of  about  1500  Hz.  It  has  been  suggested  that  grav¬ 
itational  wave  radiation  during  the  accretion  phase 
halts  the  spin  period  above  a  certain  level  (30, 31). 
The  RLDP  effect  presented  here  is  a  promising 
candidate  for  an  alternative  mechanism,  in  case  a 
sub-millisecond  AXMSP  is  detected  (8). 
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Revealing  the  Superfluid  Lambda 
Transition  in  the  Universal 
Thermodynamics  of  a  Unitary  Fermi  Gas 

Mark  ].  H.  Ku,  Ariel  T.  Sommer,  Lawrence  W.  Cheuk,  Martin  W.  Zwierlein* 

Fermi  gases,  collections  of  fermions  such  as  neutrons  and  electrons,  are  found  throughout  nature,  from 
solids  to  neutron  stars.  Interacting  Fermi  gases  can  form  a  superfluid  or,  for  charged  fermions,  a 
superconductor.  We  have  observed  the  superfluid  phase  transition  in  a  strongly  interacting  Fermi  gas  by 
high-precision  measurements  of  the  local  compressibility,  density,  and  pressure.  Our  data  completely 
determine  the  universal  thermodynamics  of  these  gases  without  any  fit  or  external  thermometer.  The 
onset  of  superfluidity  is  observed  in  the  compressibility,  the  chemical  potential,  the  entropy,  and  the  heat 
capacity,  which  displays  a  characteristic  lambda-like  feature  at  the  critical  temperature  TJTF  =  0.167(13). 
The  ground-state  energy  is  |  Ef  with  \  =  0.376(4).  Our  measurements  provide  a  benchmark  for 
many-body  theories  of  strongly  interacting  fermions. 


Phase  transitions  are  ubiquitous  in  nature: 
Water  freezes  into  ice,  electron  spins  sud¬ 
denly  align  as  materials  turn  into  magnets, 
and  metals  become  superconducting.  Near  the 


transitions,  many  systems  exhibit  critical  behav¬ 
ior,  reflected  by  singularities  in  thennodynamic 
quantities:  The  magnetic  susceptibility  diverges 
at  a  ferromagnetic  transition,  and  the  specific  heat 
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shows  a  jump  at  superconducting  and  superfluid 
transitions  ( 1 ,  2),  resolved  as  the  famous  lambda 
peak  in  4He  (3).  A  novel  fonn  of  superfluidity  has 
been  realized  in  trapped,  ultracold  atomic  gases  of 
strongly  interacting  fennions,  particles  with  half¬ 
integer  spin  (4-7).  Thanks  to  an  exquisite  control 
over  relevant  system  parameters,  these  gases  have 
recently  emerged  as  a  versatile  system  well  suited 
to  solve  open  problems  in  many-body  physics  (7). 

Initial  measurements  on  the  thennodynamics 
of  strongly  interacting  Fermi  gases  have  focused 
on  trap-averaged  quantities  (8-10)  in  which  the 
superfluid  transition  is  inherently  difficult  to  ob¬ 
serve.  The  emergence  of  the  condensate  of  femiion 
pairs  in  a  spin-balanced  Fenni  gas  is  accompa¬ 
nied  by  only  minute  changes  in  the  gas  density 
(J).  Quantities  that  involve  integration  of  the  den¬ 
sity  over  the  local  potential,  such  as  the  energy 
E  (11)  and  the  pressure  P  (12),  are  only  weakly 
sensitive  to  the  sudden  variations  in  the  thermo¬ 
dynamics  of  the  gas  expected  near  the  superfluid 
phase  transition  (13). 

For  a  neutral  gas,  thermodynamic  quantities 
involving  the  second  derivative  of  the  pressure  P 
are  expected  to  become  singular  at  the  second- 
order  phase  transition  into  the  superfluid  state. 
An  example  is  the  isothermal  compressibility 
K  =  nlh  ly,  the  relative  change  of  the  gas  density 
n  due  to  a  change  in  the  pressure  P.  Because  the 
change  in  pressure  is  related  to  the  change  in 
chemical  potential  p  of  the  gas  via  dP  =  n  dp  at 
constant  temperature,  k  =  T  |^  |  r  is  a  second  de¬ 
rivative  of  the  pressure,  and  thus  should  reveal 
a  clear  signature  of  the  transition. 

The  general  strategy  to  detennine  the  ther¬ 
modynamic  properties  of  a  given  substance  is  to 
measure  an  equation  of  state  (EoS),  such  as  the 
pressure  P(p,T)  as  a  function  of  the  chemical  po¬ 
tential  p  and  the  temperature  T.  Equivalently,  re¬ 
placing  the  pressure  by  the  density  n  =  @  \T,  one 
can  detennine  the  density  EoS  n(  p,T).  We  directly 
measure  the  local  gas  density  n(V)  as  a  function  of 
the  local  potential  V  from  in  situ  absorption 
images  of  a  trapped,  strongly  interacting  Fermi 
gas  of  6Li  atoms  at  a  Feshbach  resonance  (5). 
The  trapping  potential  is  cylindrically  symmetric, 
with  harmonic  confinement  along  the  axial  direc¬ 
tion;  this  symmetry  allows  us  to  find  the  three- 
dimensional  (3D)  density  through  the  inverse 
Abel  transform  of  the  measured  column  density 
(14, 15).  The  local  potential  is  directly  determined 
from  the  atomic  density  distribution  and  the  ac¬ 
curately  known  hamionic  potential  along  the  axial 
direction. 

The  compressibility  k  follows  as  the  change 
of  the  density  n  with  respect  to  the  local  potential 
^experienced  by  the  trapped  gas.  The  change  in 
the  local  chemical  potential  is  given  by  the  neg¬ 
ative  change  in  the  local  potential,  dp  =  -dV,  and 
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hence  the  local  compressibility  is  k  =  |r. 

We  can  then  replace  the  unknown  chemical  po¬ 
tential  p  in  the  density  EoS  n(\i,T)  by  the  known 
variation  of  n  with  p  in  the  atom  trap,  given  by  k. 
Instead  of  the  a  priori  unknown  temperature  T, 
we  determine  the  pressure  P(V)  =  J_^dp'n(p')  = 
\vdV'n(V)  given  by  the  integral  of  the  density 
over  the  potential  (16).  The  resulting  equation  of 
state  n( k,P)  contains  only  quantities  that  can  be 
directly  obtained  from  the  density  distribution. 
This  represents  a  crucial  advance  over  previous 
methods  that  require  the  input  of  additional  ther¬ 
modynamic  quantities,  such  as  the  temperature  T 
and  the  chemical  potential  p,  whose  determina¬ 
tion  requires  the  use  of  a  fitting  procedure  or  an 
external  thermometer,  as  in  (11,  12). 

We  normalize  the  compressibility  and  the  pres¬ 
sure  by  the  respective  quantities  at  the  same 
local  density  for  a  noninteracting  Fenni  gas  at 
T  =  0,  Kp  =  2 hEp  ,  and  P0  =  \nEp,  where 
Ef  =  h  /  is  the  local  Fermi  energy  and  m  is 
the  particle  mass,  yielding  k  =  k/k0  and p  =  P/Po- 
For  dilute  gases  at  the  Feshbach  resonance,  the 
scattering  length  diverges  and  is  no  longer  a  rele¬ 
vant  length  scale.  In  the  absence  of  an  interaction- 
dependent  length  scale,  the  thennodynamics  of 
such  resonant  gases  are  universal  (17),  and  k  is  a 
universal  function  of  p  only.  Every  experimental 
profile  n(V),  irrespective  of  the  trapping  potential, 
the  total  number  of  atoms,  or  the  temperature,  must 
produce  the  same  universal  curve  k  versus  p.  By 
averaging  many  profiles,  one  obtains  a  low-noise 
detennination  of  k  ( p). 


Our  method  has  been  tested  on  the  nonin¬ 
teracting  Fenni  gas  that  can  be  studied  in  two  in¬ 
dependent  ways:  in  spin-balanced  gases  near  the 
zero-crossing  of  the  scattering  length  and  in  the 
wings  of  highly  imbalanced  clouds  at  unitarity, 
where  only  one  spin  state  is  present  locally.  Both 
determinations  yield  the  same  noninteracting  com¬ 
pressibility  EoS  (Fig.  1). 

Figure  1  also  shows  the  compressibility  equa¬ 
tion  of  state  k  (p)  for  the  unitary  Fermi  gas.  In  the 
high-temperature  (p  »  1)  regime,  the  pressure, 
and  hence  all  other  thennodynamic  quantities,  al¬ 
low  for  a  Virial  expansion  in  terms  of  the  fugacity 
e ^  (18):  P  fWi3  =  2 JT  bj  eyPt\  with  the  nth-order 
Virial  coefficients  bn.  It  is  known  that  bi  =  1,  b2  = 
3x/2/8,  and  b3  =  -0.29095295  (18):  our  data  show 
good  agreement  with  the  third-order  Virial  expan¬ 
sion.  Fixing  b2  and  b3,  our  measurement  yields  a 
prediction  for  h4  =  +0.065(  1 0),  in  agreement  with 
(12),  but  contradicting  a  recent  four-body  calcu¬ 
lation  that  gives  a  negative  sign  (19). 

At  degenerate  temperatures  (p  <  1),  the  nor¬ 
malized  compressibility  rises  beyond  that  of  a 
noninteracting  Fermi  gas,  as  expected  for  an  at¬ 
tractively  interacting  gas.  A  sudden  rise  of  the  com¬ 
pressibility  at  around  p  =  0.55,  followed  by  a 
decrease  at  lower  temperatures  marks  the  super¬ 
fluid  transition.  The  expected  singularity  of  the 
compressibility  at  the  transition  is  rounded  off  by 
the  finite  resolution  of  our  imaging  system.  Be¬ 
low  the  transition  point,  the  decrease  of  the  com¬ 
pressibility  is  consistent  with  the  expectation  from 
Bardeen-Cooper-Schrieffer  (BCS)  theory,  in  which 


Fig.  1.  Normalized  compressibility  k/k0  versus  normalized  pressure  PIP0  of  the  unitary  Fermi  gas  (red 
solid  circles).  Each  data  point  is  the  average  of  between  30  and  150  profiles.  The  error  bars  show  mean  + 
SD,  including  systematic  errors  from  image  calibration  (13).  Blue  solid  line:  third-order  Virial  expansion. 
Black  open  squares  (black  open  diamonds):  data  for  a  noninteracting  Fermi  gas  obtained  with  a  highly  spin- 
imbalanced  mixture  at  the  Feshbach  resonance  (spin-balanced  gas  near  zero-crossing  of  the  scattering 
length).  Black  solid  curve:  theory  for  a  noninteracting  Fermi  gas.  Black  dashed  curve:  the  relation  k  =  lip 
that  must  be  obeyed  at  zero  temperature  both  for  the  noninteracting  gas  (k  =  1/p  =  1)  and  the  unitary  gas 
(k  =  1/p  =  lie)  (dotted  lines).  Gray  band:  the  uncertainty  region  for  the  T=  0  value  of  k  =  and  p  =  i;. 
Blue  dashed  curve:  model  for  the  EoS  of  the  unitary  Fermi  gas  [above  Tc:  interpolation  from  the  Monte 
Carlo  calculation  (34):  below  Tc:  BCS  theory,  including  phonon  and  pair-breaking  excitations].  Green  solid 
curve:  effect  of  2  pm  optical  resolution  on  the  model  EoS. 
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single-particle  excitations  freeze  out  and  pairs  form 
(see  model  in  Fig.  1). 

As  T  — >  0,  the  Fermi  energy  EF  is  the  only 
intensive  energy  scale,  so  the  chemical  potential 
must  be  related  to  EF  by  a  universal  number,  p  = 
\EF,  where  \  is  known  as  the  Bertsch  parameter 
( 6 ,  7).  It  follows  that  at  T=  0,  k=  Vp  =  V%  (13). 
The  extrapolation  of  the  low-temperature  exper¬ 
imental  data  for  k  (p)  toward  the  curve  k  =  1/p 
gives  i;  =  0.37(1),  a  value  that  we  find  consistently 
for  the  normalized  chemical  potential,  energy,  and 
free  energy  at  our  lowest  temperatures. 

From  the  universal  function  k  ( p ),  we  obtain 
all  other  thennodynamic  quantities  of  the  unitary 
gas.  First,  to  find  the  normalized  temperature  T/TF 
(where  kB  TF  =  EF),  note  that  the  change  in  pres¬ 
sure  with  TITp  at  constant  temperature  is  re¬ 
lated  to  the  compressibility.  One  finds  = 

(p  -  i),  so  by  integration  (13) 


where  (T/TF)f  is  the  nonnalized  temperature  at  an 
initial  normalized  pressure  p,  that  can  be  chosen 
to  lie  in  the  Virial  regime  validated  above. 


Thanks  to  the  relation  E  =  |FV,  valid  at 
unitarity  (17),  we  can  also  directly  obtain  the  heat 
capacity  per  particle  at  constant  volume  V  (13), 


CV=\8E  _  3  dp 
kgN  ~  kgNdf'N’v  ~  5d(r/7» 


Figure  2  shows  the  normalized  compressibility 
and  the  specific  heat  as  a  function  of  TITF.  At 
high  temperatures,  the  specific  heat  approaches 
that  of  a  noninteracting  Fenni  gas  and  eventually 
Cy=  |  N  kB,  the  value  for  a  Boltzmann  gas.  A 
dramatic  rise  is  observed  for  T/TF  at  around  0. 16, 
followed  by  a  steep  drop  at  lower  temperatures. 
Such  a  3,-shaped  feature  in  the  specific  heat  is 
characteristic  of  second-order  phase  transitions, 
as  in  the  famous  X.  transition  in  4He  (3).  Jumps  in 
the  specific  heat  are  well  known  from  supercon¬ 
ductors  (1)  and  3He  (2).  In  experiments  on  atomic 
gases,  such  jumps  had  only  been  inferred  from 
derivatives  to  fit  functions  that  implied  a  jump 
(20,  21).  We  do  not  expect  to  resolve  the  critical 
behavior  very  close  to  Tc.  Because  of  the  spatially 
varying  chemical  potential  in  our  trapped  sample, 
the  critical  region  is  confined  to  a  narrow  shell. 


0.2  0.4  0.6  0.8  1.0  1.2 


t/tf 


Fig.  2.  (A)  Normalized  compressibility  ic  =  [2/3]k  r?  EF  and  (B)  specific  heat  per  particle  CvINkB  of  a 
unitary  Fermi  gas  as  a  function  of  reduced  temperature  777>  (solid  red  circles).  Black  solid  curve:  theory  for 
a  noninteracting  Fermi  gas.  Blue  solid  curve:  third-order  Virial  expansion  for  the  unitary  gas.  Black  open 
squares:  data  for  the  normalized  compressibility  as  a  function  of  777>  of  a  noninteracting  Fermi  gas 
(combining  data  from  both  highly  imbalanced  gases  at  unitarity  and  balanced  gases  near  zero-crossing). 
Blue  dashed  (green  solid)  curve:  model  from  Fig.  1,  excluding  (including)  the  effect  of  finite  imaging 
resolution.  (C)  Global  condensate  fraction  at  unitarity  as  determined  from  a  rapid  ramp  to  the  molecular 
side  of  the  Feshbach  resonance,  plotted  as  a  function  of  local  777>  at  the  trap  center.  The  onset  of  con¬ 
densation  coincides  with  the  sudden  rise  of  the  specific  heat.  Error  bars,  mean  ±  SD. 
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Based  on  the  estimate  in  (22),  the  thickness  of  the 
critical  shell  is  1%  of  the  cloud  size.  The  finite 
resolution  of  our  imaging  system  (2  pm  or  about 
5%  of  the  cloud  size  in  the  radial  direction)  suf¬ 
fices  to  explain  the  rounding  of  the  singularity 
expected  from  criticality.  The  rounding  also  re¬ 
duces  the  observed  jump  in  the  heat  capacity  at 
the  transition.  We  obtain  a  lower  bound  A C/Cn  = 
(Cs  -  C„)/Cn  >  1 ,0+4,  where  CJN  (CJN)  is  the 
specific  heat  per  particle  at  the  peak  (the  onset  of 
the  sudden  rise).  Considering  the  strong  inter¬ 
actions,  this  is  surprisingly  close  to  the  BCS 
value  of  1.43  (1).  Below  Tc,  the  specific  heat  is 
expected  to  decrease  as  ~  exp(-Ao lkB  T)  due  to 
the  pairing  gap  Ao.  At  low  temperatures,  T  «TC, 
the  phonon  contribution  °c  T3  dominates  (23). 
This  behavior  is  consistent  with  our  data,  but  the 
phonon  regime  is  not  resolved. 

To  validate  our  in  situ  measurements  of  the 
superfluid  phase  transition,  we  have  employed 
the  rapid  ramp  method  to  detect  fermion  pair  con¬ 
densation  (24,  25).  The  results  (Fig.  2C)  show 
that  the  onset  of  condensation  and  the  sudden  rise 
in  specific  heat  and  compressibility  all  occur  at 
the  same  critical  temperature,  within  the  error  bars. 
Unlike  previous  experimental  determinations  of 
Tc/Tf  for  the  homogeneous  unitary  Fenni  gas 
(11, 12),  we  determine  TC/TF  directly  from  the  den¬ 
sity  profiles,  finding  a  sudden  rise  in  the  specific 
heat  and  the  onset  of  condensation  at  TC/TF  = 
0.167(13).  This  value  is  detennined  as  the  mid¬ 
point  of  the  sudden  rise,  and  the  error  is  assessed 
as  the  shift  due  to  the  uncertainty  of  the  Feshbach 
resonance  (13).  This  is  in  very  good  agreement 
with  theoretical  detenninations,  such  as  the  self- 
consistent  T-matrix  approach  that  gives  TC/TF  ~ 
0.16  (23),  and  Monte  Carlo  calculations  that  give 
Tc/Tf=  0.171(5)  (26)  and  0.152(7)  (27).  There  is 
a  current  debate  on  the  possibility  of  a  pseudo¬ 
gap  phase  of  prefonned  pairs  above  Tc  (12,  28). 
A  pairing  gap  for  single-particle  excitations  above 
the  transition  should  be  signaled  by  a  downturn 
of  the  specific  heat  above  Tc,  which  is  not  ob¬ 
served  in  our  measurements. 

From  the  definition  of  the  compressibility 
K  =  \t  > we  can  obtain  the  reduced  chemical 
potential  \i/Ef  as  a  function  of  the  TITF  (Fig.  3A) 
(13).  This  function  is  here  obtained  from  measured 
quantities,  rather  than  from  numerical  derivatives 
of  data  that  involved  uncontrolled  thermometry  (11). 
In  the  interval  of  T/TF  from  around  0.25  to  1, 
the  chemical  potential  is  close  to  that  of  a  non¬ 
interacting  Fenni  gas,  shifted  by  (£,„  -  1  )EF  be¬ 
cause  of  interactions  present  in  the  nonnal  state, 
with  4 „  ~  0.45.  Unlike  a  nonnal  Fenni  gas,  the 
chemical  potential  attains  a  maximum  of  | i/EF  = 
0.42(1)  at  T/Tf~  0.171(10),  and  then  decreases  at 
lower  temperatures,  as  expected  for  a  superfluid 
of  paired  fermions  (23).  As  the  temperature  is  in¬ 
creased  from  zero  in  a  superfluid,  first  the  emer¬ 
gence  of  phonons  (sound  excitations)  and  then 
the  breaking  of  fermion  pairs  contribute  to  in¬ 
creasing  the  chemical  potential.  At  Tc,  the  sin¬ 
gular  compressibility  implies  a  sharp  change  in 
slope  for  p !Ef,  in  agreement  with  our  observa- 
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tion  and  theory  (23).  At  low  temperatures,  the 
reduced  chemical  potential  \i/Ef  saturates  to  the 
universal  value As  the  internal  energy  E  and 
the  free  energy  F  satisfy  E(T)>  E( 0)  =| N^EF= 
F( 0)  >  F(T)  for  all  T,  the  reduced  quantities 

fE=lm^  =  P  and/^fj^  =  (Fig. 

3  A)  provide  upper  and  lower  bounds  for  \  (29). 
Taking  the  coldest  points  of  these  three  curves  and 
including  the  systematic  error  due  to  the  effective 
interaction  range,  we  find  i;  =  0.376(4).  The  un¬ 
certainty  in  the  Feshbach  resonance  is  expected 
to  shift  '%  by  at  most  2%  (13).  This  value  is  con¬ 
sistent  with  a  recent  upper  bound  ^  <  0.383(1)  from 
(30),  is  close  to  i;  =  0.36(1)  from  a  self-consistent 
T-matrix  calculation  (23),  and  agrees  with  i;  = 
0.367(9)  from  an  epsilon  expansion  (31).  It  lies 
below  earlier  estimates  E,  =  0.44(2)  (32)  and  E,  = 
0.42(1)  (33)  from  fixed-node  quantum  Monte  Carlo 
calculation  that  provides  upper  bounds  on  E,.  Our 
measurement  agrees  with  several  less  accurate  ex¬ 
perimental  determinations  (6)  but  disagrees  with 
the  most  recent  experimental  value  0.415(10)  that 
was  used  to  calibrate  the  pressure  in  (12). 

From  the  energy,  pressure,  and  chemical  po¬ 
tential,  we  can  obtain  the  entropy  S  =  j(E  +  PV  - 
( jlN),  and  hence  the  entropy  per  particle  S  /NkB  — 


II 

T 


as  a  function  of  777^.  (Fig.  3B).  At 


high  temperatures,  S  is  close  to  the  entropy  of 
an  ideal  Fermi  gas  at  the  same  TITp.  Above  Tc, 
the  entropy  per  particle  is  nowhere  small  com¬ 
pared  with  kB.  Also,  the  specific  heat  Cv  is  not 
linear  in  T  in  the  normal  phase.  This  shows  that 
the  normal  regime  above  Tc  cannot  be  described  in 
tenns  of  a  Landau  Fermi  Liquid  picture,  although 
some  thermodynamic  quantities  agree  surpris¬ 
ingly  well  with  the  expectation  for  a  Fenni  liquid 
[see  (12)  and  (73)].  Below  about  T/TF=  0.17,  the 
entropy  starts  to  strongly  fall  off  compared  with  that 
of  a  noninteracting  Fermi  gas,  which  we  again 
interpret  as  the  freezing  out  of  single-particle  excita¬ 
tions  as  a  result  of  the  formation  of  fennion  pairs. 
Far  below  Tc,  phonons  dominate.  They  only  have  a 
minute  contribution  to  the  entropy  (23),  less  than 
0.02  kB  at  TITp  =0.1,  consistent  with  our  measure¬ 
ments.  At  die  critical  point,  we  obtain  Sc  =  0.73(13) 
NkB,  in  agreement  with  theory  (23).  It  is  encourag¬ 
ing  for  future  experiments  with  fennions  in  optical 
lattices  that  we  obtain  entropies  less  than  0.04  N 
kB,  far  below  critical  entropies  required  to  reach 
magnetically  ordered  phases. 

From  the  chemical  potential  [i/EF  and  T/TF  = 
(3h2)2/3  (nx3)2/3, we  finally  obtain  the  density  EoS 

/i(p,r)=^_4,([3p),  with  the  de  Broglie  wave¬ 
length  X  =  The  pressure  EoS  follows 


as  P(\x,T)  =  ^t./p(Ph),  with  fP  =  f^p/TPp). 
Figure  4  shows  the  density  and  pressure  nor¬ 
malized  by  their  noninteracting  counterparts  at 
the  same  chemical  potential  and  temperature.  For 
the  nonnal  state,  a  concurrent  theoretical  calcu¬ 
lation  employing  a  new  Monte  Carlo  method 
agrees  excellently  with  our  data  (34).  Our  data 


deviate  from  a  previous  experimental  determi-  E,  =  0.415(10)  (35)  and  disagree  with  the  energy 
nation  of  the  pressure  EoS  (12)  that  was  cal-  measurement  in  (11)  that  used  a  thermometry  in- 
ibrated  with  an  independently  measured  value  of  consistent  with  the  Virial  expansion  (10).  Around 


Fig.  3.  (A)  Chemical  potential  p,  energy  E,  and  free  energy  F  of  the  unitary  Fermi  gas  versus  777>.  p  (red 
solid  circles)  is  normalized  by  the  Fermi  energy  EF,  and  E  (black  solid  circle)  and  F  (green  solid  circle)  are 
normalized  by  E0  =  § N  EF.  At  high  temperatures,  all  quantities  approximately  track  those  for  a  non¬ 
interacting  Fermi  gas,  shifted  by  -  1  (dashed  curves).  The  peak  in  the  chemical  potential  signals  the 
onset  of  superfluidity.  In  the  deeply  superfluid  regime  at  low  temperatures,  p/Ff,  EIE0,  and  FIF0  all  approach 
%  (blue  dashed  line).  (B)  Entropy  per  particle.  At  high  temperatures,  the  entropy  closely  tracks  that  of  a 
noninteracting  Fermi  gas  (black  solid  curve).  The  open  squares  are  from  the  self-consistent  T-matrix 
calculation  (23).  A  few  representative  error  bars  are  shown,  representing  mean  +  SD. 


Fig.  4.  (A)  Density  and  (B)  pressure  of  a  unitary  Fermi  gas  versus  p/kB  T,  normalized  by  the  density  and 
pressure  of  a  noninteracting  Fermi  gas  at  the  same  chemical  potential  p  and  temperature  T.  Red  solid 
circles:  experimental  EoS.  Blue  dashed  curves:  low-temperature  behavior  with  \  =  0.364  (upper),  0.376 
(middle),  and  0.388  (lower).  Black  dashed  curve:  low-temperature  behavior  with  \  at  upper  bound  of  0.383 
from  (30).  Green  solid  circles  (black  fine  dashed  line):  MIT  experimental  data  (theory)  for  the  ideal  Fermi 
gas.  Blue  solid  squares  (blue  curve):  diagrammatic  Monte  Carlo  calculation  (34)  for  density  (pressure,  with 
blue  dashed  curves  denoting  the  uncertainty  bands).  Solid  green  line:  third-order  Virial  expansion.  Open 
black  squares:  self-consistent  T-matrix  calculation  (23).  Open  green  circles:  lattice  calculation  (36).  Orange 
star  and  blue  triangle:  critical  point  from  the  Monte  Carlo  calculations  (26)  and  (27),  respectively.  Solid 
diamonds:  Ecole  Normale  Superieure  experiment  (12).  Purple  open  diamonds:  Tokyo  experiment  (11). 
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the  critical  point,  the  density  shows  a  strong  var¬ 
iation,  whereas  the  pressure,  the  integral  of  the 
density  over  p  at  constant  T,  is  naturally  less  sen¬ 
sitive  to  the  superfluid  transition. 

In  conclusion,  we  have  performed  thennody- 
namic  measurements  of  the  unitary  Fermi  gas 
across  the  superfluid  phase  transition  at  the  level 
of  uncertainty  of  a  few  percent,  without  any  fits 
or  input  from  theory,  enabling  validation  of  the¬ 
ories  for  strongly  interacting  matter.  Similar  un¬ 
biased  methods  can  be  applied  to  other  systems, 
for  example,  two-dimensional  Bose  and  Fermi 
gases  or  fermions  in  optical  lattices. 
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Iron  Catalysts  for  Selective 
Anti-Markovnikov  Alkene 
Hydrosilylation  Using  Tertiary  Silanes 

Aaron  M.  Tondreau,1  Crisita  Carmen  Hojilla  Atienza,1  Keith  ].  Weller,2  Susan  A.  Nye,2 
Kenrick  M.  Lewis,3  Johannes  G.  P.  Delis,4  Paul  J.  Chirik1* * 

Alkene  hydrosilylation,  the  addition  of  a  silicon  hydride  (Si-H)  across  a  carbon-carbon  double  bond, 
is  one  of  the  largest-scale  industrial  applications  of  homogeneous  catalysis  and  is  used  in  the 
commercial  production  of  numerous  consumer  goods.  For  decades,  precious  metals,  principally 
compounds  of  platinum  and  rhodium,  have  been  used  as  catalysts  for  this  reaction  class.  Despite 
their  widespread  application,  limitations  such  as  high  and  volatile  catalyst  costs  and  competing 
side  reactions  have  persisted.  Here,  we  report  that  well-characterized  molecular  iron  coordination 
compounds  promote  the  selective  anti-Markovnikov  addition  of  sterically  hindered,  tertiary  silanes 
to  alkenes  under  mild  conditions.  These  Earth-abundant  base-metal  catalysts,  coordinated  by 
optimized  bis(imino)pyridine  ligands,  show  promise  for  industrial  application. 


Metal-catalyzed  olefin  hydrosilylation,  which 
fonns  alkylsilanes  by  cleaving  a  silicon- 
hydrogen  bond  and  adding  the  frag¬ 
ments  across  a  carbon-carbon  double  bond  (1,2), 
finds  widespread  application  in  the  commercial 
manufacture  of  silicone-based  surfactants,  fluids, 
molding  products,  release  coatings,  and  pressure- 
sensitive  adhesives  (3,  4).  Consequently,  hydro¬ 
silylation  has  emerged  as  one  of  the  largest-scale 
applications  of  homogeneous  catalysis  (5-9). 


For  more  than  three  decades,  precious  metal  com¬ 
pounds  with  Pt,  Pd,  Ru,  and  Rh  have  been  used 
almost  exclusively  as  catalysts.  Platinum  com¬ 
pounds  such  as  Karstedt’s  and  Speier’s  cata¬ 
lysts,  Pt2{[(CH2=CH)SiMe2]20}3  (Me,  methyl)  and 
H2PtCl6-6H20/'PrOH('Pr,  isopropyl),  respective¬ 
ly,  are  the  most  widely  used  industrial  catalysts 
(1,  10-12),  though  they  suffer  from  chemical  lim¬ 
itations  such  as  intolerance  to  amino-substituted 
olefins  and  a  tendency  to  catalyze  competing  isom¬ 


erization  of  the  tenninal  alkenes  to  internal  isomers. 
Undesired  isomerization  often  necessitates  sub¬ 
sequent  purification  steps  that  are  both  energy 
and  cost  intensive.  Furthennore,  decomposition 
of  the  catalyst  to  colloidal  platinum  contributes  to 
unwanted  side  reactions  and  also  causes  dis¬ 
coloration  of  the  final  products. 

It  has  been  estimated  that  the  worldwide  sil¬ 
icone  industry  consumed  -180,000  troy  ounces 
(5.6  metric  tons)  of  platinum  in  2007  and  most  is  not 
recovered  (13).  The  high  cost,  coupled  with  the  in¬ 
creasing  demands  on  precious  metals  due  to  fuel-cell 
and  other  emerging  technologies,  has  increased  the 
volatility  of  the  platinum  market  (14).  The  combi¬ 
nation  of  chemical,  economic,  and  political  chal¬ 
lenges  inspires  the  exploration  of  inexpensive  and 
Earth-abundant  catalysts  using  iron,  manganese, 
and  cobalt  (15).  At  the  core  of  this  challenge  is  sup¬ 
pressing  tendencies  of  first-row  transition  metals 
toward  one-electron  redox  processes  in  favor  of  the 
two-electron  chemistry  associated  with  the  heavier 
metals  that  probably  make  up  the  fundamental 
steps  in  a  catalytic  cycle  for  alkene  hydrosilylation. 
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Quantum  degenerate  Bose-Fermi  mixture  of  chemically  different  atomic  species 

with  widely  tunable  interactions 

Jee  Woo  Park,1,2  Cheng-Hsun  Wu,1'2  Ibon  Santiago,1,2  Tobias  G.  Tiecke,1,3  Sebastian  Will,1,2 
Peyman  Ahmadi,1-2  and  Martin  W.  Zwierlein1-2 
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We  have  created  a  quantum  degenerate  Bose-Fermi  mixture  of  23Na  and  40K  with  widely  tunable  interactions 
via  broad  interspecies  Feshbach  resonances.  Over  30  Feshbach  resonances  between  23Na  and  40K  were  identified, 
including  p-wave  multiplet  resonances.  The  large  and  negative  triplet  background  scattering  length  between  23Na 
and  40K  causes  a  sharp  enhancement  of  the  fermion  density  in  the  presence  of  a  Bose  condensate.  As  explained  via 
the  asymptotic  bound-state  model,  this  strong  background  scattering  leads  to  wide  Feshbach  resonances  observed 
at  low  magnetic  fields.  Our  work  opens  up  the  prospect  to  create  chemically  stable,  fermionic  ground-state 
molecules  of  23Na-40K,  where  strong,  long-range  dipolar  interactions  would  set  the  dominant  energy  scale. 

DOI:  10. 1 103/PhysRevA.85.05 1602  PACS  number(s):  67.85.-d,  03.75. Ss,  34.50.-s,  37.10.De 


Ultracold  quantum  gases  realize  paradigms  of  condensed- 
matter  physics  in  pristine  fashion,  such  as  the  superfluid 
to  Mott  insulator  transition  [1],  the  BEC-BCS  crossover  in 
fermionic  superfluids  [2,3],  and  the  Berezinskii-Kosterlitz- 
Thouless  transition  in  two-dimensional  Bose  gases  [4],  A 
plethora  of  novel  many-body  systems  may  become  accessible 
through  the  advent  of  quantum  mixtures  of  different  atomic 
species.  In  particular,  Bose-Fermi  mixtures  with  widely 
tunable  interactions  should  reveal  boson-mediated  interactions 
between  fermions  and  possibly  boson-induced  p- wave  super¬ 
fluidity  [5,6].  The  fate  of  impurities  in  a  Fermi  sea  [7]  or  a  Bose 
condensate  [8-10]  can  be  studied,  and  new  quantum  phases  of 
matter  are  predicted  in  optical  lattices  [11].  Furthermore,  the 
creation  of  fermionic  ground- state  molecules  starting  from  a 
degenerate  Bose-Fermi  mixture  opens  up  an  intriguing  avenue 
of  research,  as  this  results  in  a  Fermi  gas  with  long-range, 
anisotropic  dipole-dipole  interactions  [12],  Since  the  first 
degenerate  Bose-Fermi  mixture  of  different  atomic  species, 
23Na  and  6Li  [13],  a  variety  of  such  systems  has  been  realized 
[8,14-23],  However,  so  far  only  one  mixture,  87Rb-40K,  has 
allowed  tunability  of  interspecies  interactions  with  relative 
ease  by  means  of  a  moderately  wide  (A  B  ~  3  G)  Feshbach 
resonance  [24],  and  only  in  this  case  fermionic  Feshbach 
molecules  have  successfully  been  produced  [25,26]. 

In  this  Rapid  Communication,  we  report  on  the  experimen¬ 
tal  realization  of  a  unique  Bose-Fermi  mixture  of  23Na  and 
40K  and  the  observation  of  over  30  s-  and  p- wave  Feshbach 
resonances  at  low  magnetic  fields.  We  demonstrate  that  23Na 
is  an  efficient  coolant  for  sympathetic  cooling  of  40K.  A 
pattern  of  wide  .v-wave  resonances  exists  for  most  of  the 
energetically  stable  hyperfine  combinations,  the  widest  being 
located  at  138  G  with  a  width  of  about  30  Gin  the  23Na|F  =  1, 
mp  =  1)  +40K|U  =  9/2 ,mp  =  —5/2)  hyperfine  configura¬ 
tion.  We  also  observe  p-wave  multiplet  resonances  that 
are  resolved  thanks  to  their  location  at  low  magnetic 
fields. 

In  the  singlet  rovibrational  ground  state,  the  NaK  molecule 
is  known  to  have  a  large  permanent  electric  dipole  moment  of 


2.72(6)  D  [27,28] ,  five  times  larger  than  that  of  KRb  [12],  and  is 
predicted  to  be  chemically  stable  against  atom-atom  exchange 
reactions  [29],  in  contrast  to  KRb  [30].  An  ultracold  gas  of 
fermionic  ground-state  molecules  of  NaK  will  thus  be  an  ideal 
system  for  the  study  of  Fermi  gases  with  strong,  long-range 
dipolar  interactions.  Indeed,  the  interaction  energy  here  can  be 
expected  to  be  on  the  order  of  the  Fermi  energy. 

The  experimental  setup  is  based  on  the  apparatus  presented 
in  Ref.  [8],  which  employs  two  independent  Zeeman  slowers 
capable  of  simultaneously  loading  sodium  and  potassium 
atoms  directly  into  a  UHV  chamber.  The  potassium  slower 


FIG.  1.  (Color  online)  Simultaneous  quantum  degeneracy  of  23Na 
and  40K  atoms,  (a)  and  (b),  and  (d)  and  (e)  are  pairs  of  time-of-flight 
(TOF)  absorption  images  of  a  23Na  BEC  and  a  40K  Fermi  cloud 
with  different  atom  number  balance.  A  strong  attractive  interaction 
between  the  two  species  is  observed  in  (e)  as  a  sharp  increase  of 
the  central  density  in  the  fermionic  cloud  in  the  presence  of  a  Bose 
condensate,  (c)  and  (f)  are  the  center-sliced  column  density  of  the 
fermionic  clouds  of  (b)  and  (e),  respectively.  TOF  was  (a)  11  ms, 
(b)  7  ms,  (d)  17  ms,  and  (e)  5  ms.  Atom  numbers  are  (a)  2.4  x  105, 
(b)  2  x  105,  (d)  8.2  x  105,  and  (e)  6.7  x  104. 
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selectively  decelerates  the  40K  isotope,  loading  107  atoms  into 
a  magneto-optical  trap  (MOT)  in  10  s.  For  23Na,  we  use  a 
dark  spot  MOT  [31],  which  allows  us  to  load  approximately 
109  atoms  in  2  s. 

Multispecies  experiments  can  suffer  from  atom  losses  due 
to  light-assisted  collisions  in  the  MOT  and  spin-changing 
collisions  in  the  magnetic  trap.  In  order  to  minimize  such 
losses,  we  developed  a  shelving  technique  where  23Na  is  first 
loaded  into  the  MOT,  optically  pumped  to  the  1 2, 2)  stretched 
state,  and  captured  in  the  magnetic  trap.  Next,  the  trap  gradient 
is  reduced  to  7.7  G/cm  to  only  support  the  stretched  state 
of  23Na  against  gravity.  With  this  gradient  left  on,  and  the 
23Na  thus  “shelved  in  the  dark”  in  a  purely  magnetic  trap, 
finally  the  MOT  and  slower  beams  for  40K  are  switched 
on  to  load  the  40K  MOT.  This  scheme  guarantees  that  only 
23Na  atoms  in  the  1 2, 2)  state  are  present  in  the  magnetic 
trap  before  loading  40K,  and  it  also  potentially  reduces  light- 
assisted  collisions  that  would  be  encountered  in  a  double¬ 
species  MOT. 

Once  both  species  are  loaded  into  the  optically  plugged 
magnetic  trap  [8],  the  mixture  is  cooled  for  7  s  by  rf-induced 
evaporation  of  23Na,  where  thermally  excited  23Na  atoms  in  the 
1 2, 2)  state  are  removed  from  the  trap  by  coupling  to  the  high- 
field  seeking  state  of  1 1 , 1 ) .  We  decompress  the  initial  magnetic 
gradient  of  220  G/cm  to  27.5  G/cm  at  the  end  of  evaporation 
to  reduce  three-body  losses.  The  5-//K  cold  mixture  is  then 
loaded  into  a  crossed  optical  dipole  trap  (laser  wavelength 
1064  nm,  maximum  power  4.7  and  17  W,  waist  73  and 
135  /rm). 

At  this  stage,  the  1  /e  lifetime  of  the  mixture,  with  23Na  and 
40K  still  in  their  stretched  states,  is  about  r  =  250  ms,  already 
signaling  a  strong  attractive  interaction  increasing  three-body 
losses  and  spin-changing  dipolar  losses.  We  thus  immediately 
transfer  23Na  atoms  into  their  hyperfine  ground  state  1 1 , 1 )  via  a 
Landau-Zener  sweep,  and  remove  any  remaining  |2,2)  atoms 
via  a  resonant  light  pulse.  In  the  23Na|  1,1)  +  40K|9/2,9/2) 
state,  the  mixture  now  lives  for  r  =  20  s.  The  gas  is  further 
evaporatively  cooled  in  this  spin  mixture  for  2  s  by  reducing 
the  intensity  of  the  dipole  trap  beams. 

At  the  end  of  evaporation,  a  degenerate  Fermi  gas  of 
40K  with  2  x  105  atoms  and  T/TF  =  0.6  coexists  with  a 
Bose-Einstein  condensate  of  23Na.  Two  sets  of  absorption 
images  of  the  mixture  for  different  values  of  23Na  and  40K 
atom  numbers  are  shown  in  Fig.  1,  where  atom  numbers 
are  varied  by  changing  the  MOT  loading  times  of  the  two 
species.  The  strong  attractive  interaction  between  23 Na  and 
40K  in  the  23Na|  1, 1)  +  40K|9/2,9/2)  state  is  apparent.  As  the 
condensate  grows,  the  fermionic  cloud  acquires  a  bimodal 
density  distribution  as  it  experiences  the  strong  mean-field 
potential  of  the  bosons  [25] — see  Fig.  1(e). 

The  evolution  of  the  phase  space  densities  (PSDs)  and  atom 
numbers  N  of  both  species  during  evaporation  is  shown  in 
Fig.  2  .  Temperature  is  determined  by  fitting  a  thermal  profile 
to  the  wings  of  the  23Na  cloud.  The  cooling  efficiency  F  = 
— d  ln(PSD)/<7  ln(  V)  for  sodium  in  the  magnetic  trap  is  lAa  = 
2.7,  a  rather  high  value  thanks  to  the  steep  confinement  in  the 
plugged  trap.  This  efficiency  is  not  affected  by  the  presence 
of  the  relatively  small  admixture  of  40K.  Sympathetic  cooling 
is  less  efficient  than  in  other  mixtures  [8,32]  as  40K  is  seen  to 
be  lost  due  to  three-body  collisions  in  the  magnetic  trap.  We 


10 


10 


10  ' 
Qio'3 

CO  -4 
CL  10 


-5 


10 


10  ■ 


10  ■ 
10'8  -. 


l 


»  #■ 

ODT 


■V 


T> 


CL 


Decompresed 
Magnetic  Trap 


Magnetic  Trap 


10 


10 


10  10 

Atom  Number 


5 

G> 


FIG.  2.  (Color  online)  Evolution  of  the  phase  space  density  (PSD) 
with  atom  number  ( N ).  Blue  circles:  4nK;  red  squares:  23Na  with  40K; 
green  diamonds:  23Na  alone.  Open  and  solid  symbols  represent  the 
PSD  in  the  magnetic  trap  and  optical  trap,  respectively. 


find  Tk  =  4.6  for  40K.  In  the  crossed  optical  dipole  trap,  with 
sodium  in  1 1 , 1 ) ,  the  sodium  cooling  becomes  less  efficient  due 
to  the  weaker  confinement,  T^  =  0.9,  but  the  40K  number 
is  essentially  conserved  in  this  mixture  so  that  sympathetic 
cooling  is  highly  efficient,  with  rK  =  15.3.  The  lowest  T /TF 
achieved  for  40K  after  evaporating  all  of  23Na  is  T  /  7>  =  0.35 
with  3  x  105  atoms. 

With  this  degenerate  Bose-Fermi  mixture  at  our  disposal, 
the  natural  next  step  is  to  search  for  interspecies  Feshbach 
resonances  between  23Na  and  40K.  There  have  been  theoretical 
indications  for  resonances  below  100  G  [34].  In  addition, 
a  large  and  negative  triplet  scattering  length  was  predicted 
in  Ref.  [35],  a  value  that  has  recently  been  refined  to  a,  = 
-575+!?'  [34],  indicating  that  the  triplet  potential  has  an 
almost  bound,  virtual  state  right  above  threshold.  A  large 
background  scattering  length  is  often  a  catalyst  for  wide 
Feshbach  resonances  [36],  caused  by  strong  coupling  of 
molecular  states  to  the  almost  resonant  open  channel.  A 
famous  example  is  the  300-G-wide  Feshbach  resonance  in 
6 Li  [37], 

We  performed  Feshbach  loss  spectroscopy,  mapping  out 
atom  losses  of  both  species  as  a  function  of  magnetic  field. 
Over  30  Feshbach  resonances  were  observed  in  four  different 
spin  state  combinations  of  23 Na 1 1 , 1 )  +  40 K 1 9/2 ,  m  F ) ,  from  the 
ground  spin  state  mF  =  —9/2  up  to  mF  =  —3/2.  Spin  states 
of  40K  are  prepared  starting  from  mF  —  +9/2  by  a  single 
Landau-Zener  sweep  through  the  intermediate  mF  states  at 
15  G.  The  experimentally  observed  resonance  positions  and 
widths  are  reported  in  Table  I.  Many  wide  s-wave  Feshbach 
resonances  at  low  magnetic  fields  are  identified,  the  widest  one 
at  138  G  for  collisions  between  23Na|  1, 1)  +  40K|9/2,  -  5/2), 
with  a  width  of  30  G — see  Fig.  3. 

p- wave  Feshbach  resonances  are  known  to  split  into  a 
doublet  structure  due  to  different  projections  of  the  orbital 
angular  momentum  onto  the  magnetic  field  axis  [38].  In 
the  NaK  system,  however,  we  observe  triplet  features  for 
many  p- wave  resonances — see  Fig.  4.  These  originate  from 
the  magnetic  dipole-dipole  interaction  of  constituent  atoms, 
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TABLE  I.  Data  summary  on  the  Feshbach  resonances  between  23Na  in  1 1 , 1)  and  40K  in  |9/2,mF).  The  positions  and  widths  of  the 

,  _ /  D _ D®XPt\2  i,  n2 

resonances,  B0  and  ABexpt,  are  determined  by  phenomenological  Gaussian  fits  (oce  0  '  exP‘)  to  the  observed  loss  features.  For  p- wave 

resonances,  we  report  the  positions  and  widths  of  the  resolved  features,  i.e.,  doublets  and  multiplets.  Bgh  and  ABn  give  the  theoretical  positions 
and  widths  from  the  ABM  model. 


Collision  channel 

BoXpt  (G) 

ASexp,  (G) 

K  (G) 

Afith  (G) 

Resonance  type 

23Na|  1,1)  +40K|9/2,  -9/2) 

6.35,  6.41,6.47,  6.68 

0.02 

7.2 

P 

19.12,  19.18,  19.27 

0.02 

18.3 

P 

78.3 

1.1 

75.5 

1.1 

s 

88.2 

4.3 

84.5 

5.4 

s 

23Na|  1, 1)  +  40K|9/2,  -7/2) 

7.32,  7.54 

0.2,  0.03 

8.7 

p 

23.19,  23.29 

0.05,  0.05 

22.1 

p 

81.6 

0.2 

82.1 

0.04 

s 

89.8 

1.1 

87.3 

0.6 

s 

108.6 

6.6 

105.7 

13.1 

s 

23Na|  1, 1)  +  40K|9/2,  -5/2) 

9.23,  9.60 

0.14,  0.11 

11.0 

p 

29.19,  29.45,  29.52 

0.04 

27.8 

p 

96.5 

0.5 

97.2 

0.04 

s 

106.9 

1.8 

103.8 

0.45 

s 

148  (138a) 

37  (30a) 

137.1 

26 

s 

23Na|  1, 1)  +  40K|9/2,  -3/2) 

12.51,  12.68 

0.16,  0.06 

14.8 

p 

39.39,  39.86 

0.15,  0.14 

37.2 

p 

116.9 

0.5 

118.3 

0.07 

s 

129.5 

4.6 

127.2 

0.39 

s 

175 

20 

187.8 

50.5 

s 

“The  resonance  position  and  width  have  been  refined  by  measuring  the  molecular  binding  energies  via  rf  spectroscopy  [33], 


which  couple  molecular  states  with  different  total  internal 
spin.  The  diagonal  terms  of  the  magnetic  dipole-dipole 
interaction  induce  an  energy  shift  that  differs  for  the  /«/  =  0 
and  \m/\  =  1  quantum  numbers,  giving  rise  to  the  well-known 
doublet  structures.  The  off-diagonal  terms  in  the  dipole-dipole 
interaction  couple  different  values  of  w;  while  conserving 
the  total  angular  momentum  m/  +  Mp,  where  Mp  =  m  p ,k  + 
mF  Na-  These  terms  are  in  most  mixtures  negligible  since 
molecular  states  with  different  values  of  Mp  have  to  be  nearly 
degenerate  to  result  in  a  significant  energy  shift.  However, 
due  to  the  low-field  nature  of  the  NaK  p- wave  resonances, 
multiple  molecular  states  are  nearly  degenerate  with  the 
open-channel  spin  state,  allowing  for  the  triplet  structure  to  be 
resolved. 

The  assignment  of  x-  and  p- wave  characters  of  the 
resonances  follows  from  a  simple,  but  powerful  model  of  the 
molecular  states  involved.  The  singlet  and  triplet  potentials 
of  the  interatomic  potential  allow  for  a  variety  of  bound 
states.  From  the  known  scattering  length  [34]  and  the  van 
der  Waals  coefficient  Q,  [39],  the  weakest  .v-wave  bound 
states  are  expected  at  about  Ess  —  —150(10)  MHz  and  Est  — 
—  1625(50)  MHz  for  the  singlet  and  triplet  x-wave  potentials, 
where  the  errors  reflect  uncertainties  of  the  scattering  lengths 
from  Ref.  [34].  The  p- wave  bound  states  follow  from  the 
x-wave  bound  states  as  in  Ref.  [40]  and  are  slightly  adjusted 
to  fit  the  observed  resonance  positions.  As  described  in 
Refs.  [41,42],  as  a  first  guess  for  locations  of  Feshbach 
resonances  one  can  neglect  the  part  of  the  hyperfine  interaction 
that  couples  singlet  and  triplet  bound  states.  This  already 
provides  the  pattern  of  Feshbach  resonances  as  positions 
where  these  (uncoupled)  molecular  states  cross  the  atomic 


(open-channel)  threshold.  The  analysis  shows  that  the  ob¬ 
served  resonances  are  caused  by  the  triplet  bound  states. 
Next,  we  use  the  asymptotic  bound-state  model  (ABM)  to 
include  the  singlet-triplet  coupling  of  molecular  states  [43]. 
To  couple  the  molecular  states  to  the  open  channel,  we  follow 
Marcelis  et  al.  [36]  and  only  include  the  effect  of  the  virtual 
state  causing  the  large  and  negative  triplet  scattering  length. 
The  spin  part  of  the  coupling  matrix  element  is  obtained 
from  the  ABM  Hamiltonian  and  the  spatial  part,  i.e.,  the 
wave-function  overlap  between  the  respective  bound  state  and 
the  virtual  state,  is  taken  as  one  free  fit  parameter.  For  the 
background  scattering  length  of  the  low-field  resonances  the 
effect  of  broad  resonances  is  included.  The  virtual  state  causes 
strong  coupling  of  several  x-wave  molecular  states  to  the  open 
channel,  leading  to  wide,  open-channel-dominated  resonances 
as  known  from  the  case  of  6Li.  The  theoretical  values  obtained 
with  this  model  are  shown  in  Table  I.  An  exceptionally  broad 
resonance  for  mp  —  —3/2  is  predicted  to  be  even  wider  and 
to  be  shifted  further,  possibly  hinting  at  a  shift  between  the 
loss  maximum  and  the  actual  Feshbach  resonance  position. 
Our  approach  leads  to  a  refined  triplet  bound-state  energy  of 
E*  =  -1654(3)  MHz  and  Ef  =  -1478(7)  MHz,  and  using 
the  long-range  potential  from  Ref.  [34],  we  obtain  a  refined 
value  of  the  triplet  scattering  length  of  a,  =  —  830(70)ao- 
The  errors  correspond  to  one  standard  deviation  of  a  least- 
squares  fit  to  the  eight  narrowest  x-wave  resonances  that 
are  least  sensitive  to  the  coupling  to  the  scattering  states 
and  hence  the  error  induced  by  the  ABM  is  expected  to  be 
small  [43]. 

In  conclusion,  we  have  produced  a  degenerate  Bose-Fermi 
mixture  of  23Na  and  40K,  and  identified  over  30  x-  and 
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(C)  B  (G) 


FIG.  3.  (Color  online)  Feshbach  resonances  in  23Na-40K,  here 
for  23Na|  1, 1)  +  40K|9/2,— 5/2)  collisions,  (a)  Scattering  length  from 
the  ABM  model,  (b)  Open-channel  threshold  energy  (black  solid 
line),  uncoupled  (s-wave:  dashed  lines;  p- wave:  light  blue/light 
gray  lines)  and  coupled  molecular  states  G-wave:  red/gray  solid 
lines;  p- wave:  blue/dark  gray  solid  lines).  The  blue  dots  denote 
experimentally  measured  resonances.  The  inset  shows  the  energy, 
relative  to  threshold,  of  the  molecular  state  at  the  wide  .s-wave 
resonance  at  138  G  [33].  (c)  Experimental  loss  spectra  of  40K  in 
the  presence  of  23Na.  Three  s-wave  resonances  and  two  p- wave 
manifolds  are  found,  the  latter  resolved  in  one  doublet  and  one 
triplet. 


p- wave  interspecies  Feshbach  resonances,  including  several 
exceptionally  broad  resonances.  Remarkably,  many  p- wave 
Feshbach  resonances  are  observed  to  be  triplets  or  even 


Magnetic  Field  (G) 

FIG.  4.  (Color  online)  Triplet  structure  of  the  p- wave  resonance 
at  19.1  G  for  the  23Na|l,l)  +  40K|9/2,  —  9/2)  spin  configuration.  A 
phenomenological  triple  Gaussian  fit  is  applied  as  a  guide  to  the  eyes. 


multiplets.  Our  strongly  interacting  23Na-40K  mixture  near 
these  Feshbach  resonances  should  allow  the  study  of  Bose  or 
Fermi  polarons  [7],  of  boson-mediated  interactions  between 
fermions,  and  possibly  of  unique  states  of  matter  in  optical 
lattices.  The  formation  of  fermionic  Feshbach  molecules  is 
within  reach.  In  the  rovibrational  ground  state,  NaK  molecules 
possess  a  large  induced  electric  dipole  moment  and  are  stable 
against  exchange  reactions.  One  can  thus  hope  to  create 
a  Fermi  gas  of  polar  molecules  with  strong  dipole-dipole 
interactions  that  dominate  the  many-body  physics  of  the  gas, 
rather  than  being  a  small  perturbative  effect. 
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We  follow  the  evolution  of  fermion  pairing  in  the  dimensional  crossover  from  three-dimensional  to  two- 
dimensional  as  a  strongly  interacting  Fermi  gas  of  6Li  atoms  becomes  confined  to  a  stack  of  two- 
dimensional  layers  formed  by  a  one-dimensional  optical  lattice.  Decreasing  the  dimensionality  leads  to 
the  opening  of  a  gap  in  radio-frequency  spectra,  even  on  the  Bardeen-Cooper-Schrieffer  side  of  a 
Feshbach  resonance.  The  measured  binding  energy  of  fermion  pairs  closely  follows  the  theoretical 
two-body  binding  energy  and,  in  the  two-dimensional  limit,  the  zero-temperature  mean-field  Bose- 
Einstein-condensation  to  Bardeen-Cooper-Schrieffer  crossover  theory. 
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Interacting  fermions  in  coupled  two-dimensional  (2D) 
layers  present  unique  physical  phenomena  and  are  central 
to  the  description  of  unconventional  superconductivity  in 
high-transition-temperature  cuprates  [1]  and  layered  or¬ 
ganic  conductors  [2],  Experiments  on  ultracold  gases  of 
fermionic  atoms  have  allowed  access  to  the  crossover  from 
Bose-Einstein  condensation  (BEC)  of  tightly  bound  fer¬ 
mion  pairs  to  Bardeen-Cooper-Schrieffer  (BCS)  superflu¬ 
idity  of  long-range  Cooper  pairs  in  three  spatial 
dimensions  [3,4]  and,  more  recently,  the  confinement  of 
interacting  Fermi  gases  to  two  spatial  dimensions  [5-9],  A 
fermionic  superfluid  loaded  into  a  periodic  potential  should 
form  stacks  of  two-dimensional  superfluids  with  tunable 
interlayer  coupling  [10-13],  an  ideal  model  for  Josephson- 
coupled  quasi-2D  superconductors  [1,14].  For  deep  poten¬ 
tials  in  the  regime  of  uncoupled  2D  layers,  increasing  the 
temperature  of  the  gas  is  expected  to  destroy  superfluidity 
through  the  Berezinskii-Kosterlitz-Thouless  mechanism 
[15-17],  while  more  exotic  multiplane  vortex  loop  excita¬ 
tions  are  predicted  for  a  three-dimensional  (3D)  aniso¬ 
tropic  BCS  superfluid  near  the  critical  point  [18]. 

In  this  Letter,  we  study  fermion  pairing  across  the  cross¬ 
over  from  3D  to  2D  in  a  periodic  potential  of  increasing 
depth.  To  form  a  bound  state  in  3D,  the  attraction  between 
two  particles  in  a  vacuum  must  exceed  a  certain  threshold. 
However,  if  the  two  particles  interact  in  the  presence  of  a 
Fermi  sea,  the  Cooper  mechanism  allows  pairing  for  arbi¬ 
trarily  weak  interactions  [19].  In  2D,  even  two  particles  in  a 
vacuum  can  bind  for  arbitrarily  weak  interactions. 
Surprisingly,  the  mean-field  theory  of  the  BEC-BCS  cross¬ 
over  in  2D  predicts  that  the  binding  energy  of  fermion  pairs 
in  the  many-body  system  is  identical  to  the  two-body 
binding  energy  Eb  [20].  Indeed,  to  break  a  pair  and  remove 
one  pairing  partner  from  the  system  costs  an  energy  [21] 

£),  mf  =  V p2  +  A2  —  j±  within  mean-field  theory,  where 
/x  is  the  chemical  potential  and  A  is  the  pairing  gap.  In  2D, 
one  finds  [20]  p.  =  EF  —  Eh/2  and  A2  =  2 EFEb,  where 
Ef  is  the  Fermi  energy,  and  thus  EbMP  =  Eh ;  i.e.,  the 


PACS  numbers:  67.85. -d,  03.75.Lm,  05.30.Fk,  32.30.Bv 

many-body  and  two-body  binding  energies  are  predicted 
to  be  identical  throughout  the  BEC-BCS  crossover. 

We  realize  a  system  that  is  tunable  from  3D  to  2D  with  a 
gas  of  ultracold  fermionic  6Li  atoms  trapped  in  an  optical 
trap  and  a  standing-wave  optical  lattice.  The  lattice  pro¬ 
duces  a  periodic  potential  along  the  z  direction, 

V(z)  =  V0si  rr{Trz/d),  (1) 

with  depth  V0  and  lattice  spacing  d  =  532  nm.  Together 
with  the  optical  trap,  the  lattice  interpolates  between  the 
3D  and  2D  limits.  It  gradually  freezes  out  motion  along  one 
dimension  and  confines  particles  in  increasingly  uncoupled 
layers.  Features  characteristic  of  the  2D  system  appear  as 
the  strength  of  the  periodic  potential  is  increased.  The 
threshold  for  pairing  is  reduced,  allowing  pairs  to  form 
for  weaker  attractive  interactions  than  in  the  3D  system. 
The  effective  mass  of  particles  increases  along  the  confined 
direction,  and  the  center  of  mass  and  relative  degrees  of 
freedom  of  an  atom  pair  become  coupled  [11].  For  a  deep 
potential  that  suppresses  interlayer  tunneling,  the  system 
is  an  array  of  uncoupled  two-dimensional  layers.  Here, 
the  center  of  mass  and  relative  motion  decouple  and  fer¬ 
mion  pairs  form  for  the  weakest  interatomic  attraction 
[11,22,23], 

In  the  experiment,  the  appearance  of  bound  fermion 
pairs  is  revealed  using  radio-frequency  (rf)  spectroscopy. 
The  atomic  gas  consists  of  an  equal  mixture  of  6Li  atoms  in 
the  first  and  third  hyperfine  states  (denoted  as  |1)  and  |3)), 
chosen  to  minimize  final-state  interaction  effects  in  the  rf 
spectra  [24],  Interactions  between  atoms  in  states  |1)  and 
1 3)  are  greatly  enhanced  by  a  broad  Feshbach  resonance  at 
690.4(5)  G  [25].  An  rf  pulse  is  applied  to  transfer  atoms 
from  one  of  the  initial  hyperfine  states  to  the  unoccupied 
second  hyperfine  state  (denoted  as  |2)).  In  previous  work 
on  rf  spectroscopy  of  40K  fermions  in  a  deep  one¬ 
dimensional  (ID)  lattice  [8],  an  rf  pulse  transferred  atoms 
from  an  initially  weakly  interacting  state  into  a  strongly 
interacting  spin  state,  likely  producing  polarons  [26].  In 
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our  Letter,  the  initial  state  is  the  strongly  interacting, 
largely  paired  Fermi  gas  in  equilibrium,  and  the  final  state 
is  weakly  interacting. 

An  asymmetric  dissociation  peak  (the  bound-to-free 
transition)  in  the  rf  spectrum  indicates  the  presence  of 
fermion  pairs.  For  two-particle  binding,  the  pair  dissocia¬ 
tion  line  shape  in  the  3D  and  2D  limits  is  proportional  to 
p(hv  —  Eh)/v2,  with  p  the  free-particle  density  of  states 
and  v  =  ±{vrf  —  vM)  the  offset  of  the  rf  frequency  nrl 
from  the  hyperfine  splitting  vM  (plus  symbol:  1 1)  — »  |2) 
transition;  minus  symbol:  |3)  — > ►  |2)  transition).  This  form 
can  be  obtained  from  Fermi’s  golden  rule  and  the  bound- 
state  wave  function  in  momentum  space;  see  also 
Refs.  [21,27].  In  2D,  the  expected  dissociation  line  shape 
is  then  proportional  to 


I{v) 


0(hv  -  Eb) 


(2) 


In  addition  to  the  pairing  peak,  at  finite  temperature  one 
expects  a  peak  in  the  rf  spectrum  due  to  unbound  atoms 
(the  free-to-free  transition).  A  narrow  bound- to-bound 
transition  can  also  be  driven  at  an  offset  frequency  vhh  = 
(Eb  —  E'h)/h  that  transfers  one  spin  state  of  the  initial 
bound  pair  with  binding  energy  Eb  into  a  bound  state  of 
1 2)  with  |1)  or  |3),  of  binding  energy  E'b.  For  a  1 1)  —  |3) 
mixture  near  the  Feshbach  resonance,  Eb  <5C  E'h  [24],  so 
the  bound-to-bound  peak  is  well- separated  from  the 
bound-to-free  and  free-to-free  peaks.  As  very  recently 
calculated  [28],  final-state  interactions  and  the  anomalous 
nature  of  scattering  in  2D  introduce  an  additional  factor  of 

i„2/ r  J  pi  \ 

In  1[(hv-Eb)/E'h]+ni  int0  E T  (2)>  causinS  a  rounding  off  of  the 
sharp  peak  expected  from  the  step  function. 

In  a  ID  lattice,  the  binding  energy  for  two-body  pairs  is 
determined  by  the  lattice  spacing  d,  the  depth  Vo,  and  the 

3D  scattering  length  a.  In  the  2D  limit  V0  »  ER,  with 

*,2  2 

recoil  energy  ER  =  the  scattering  properties  of  the 
gas  are  completely  determined  by  Eb  [22,23].  In  that  limit, 
the  lattice  wells  can  be  approximated  as  harmonic  traps 
with  level  spacing  hcoz  =  2-JV(jEK  and  harmonic  oscillator 

length  lz  =  In  a  many-particle  system  in  2D,  the 

ratio  of  the  binding  energy  to  the  Fermi  energy  determines 
the  strength  of  interactions.  The  2D  scattering  amplitude 
f =  -ln(J2I)  +  ,V/2  for  collisions  with  energy  EF  is 
parametrized  by  InlA'/.noo),  where  kF  =  , J2mEF/h  and 
fl2D  =  h/y/mEb.  It  is  large  when  |  \n(kFa2o)\  £  1 
[22,23],  corresponding  to  the  strong-coupling  regime 
[28,29].  The  BEC  side  of  the  BEC-BCS  crossover  corre¬ 
sponds  to  negative  values  of  \n(kFa2F)),  while  the  BCS  side 
corresponds  to  positive  values  [20]. 

The  experimental  sequence  proceeds  as  follows.  An 
ultracold  gas  of  6Li  is  produced  by  sympathetic  cooling 
with  ;iNa  as  described  previously  [21].  The  6Li  atoms  are 
transferred  from  a  magnetic  trap  to  an  optical  dipole  trap 
(wavelength  1064  nm,  waist  120  pu m),  with  axial  harmonic 


confinement  (frequency  22.8  Hz)  provided  by  magnetic 
field  curvature.  With  6Li  polarized  in  state  |1),  the  mag¬ 
netic  bias  field  is  raised  to  568  G,  and  an  equal  mixture  of 
hyperfine  states  |1)  and  |3)  is  created  using  a  50%  rf 
transfer  from  |1)  to  |2)  followed  by  a  full  transfer  from 
1 2)  to  |3).  The  field  is  then  raised  to  the  final  value,  and 
evaporative  cooling  is  applied  by  lowering  the  depth  of  the 
optical  dipole  trap,  resulting  in  a  fermion  pair  condensate 
with  typically  5  X  105  atoms  per  spin  state.  The  lattice  is 
then  ramped  up  over  100  ms.  The  retro-reflected  lattice 
beam  (wavelength  1064  nm)  is  at  an  angle  of  0.5  degrees 
from  the  optical  dipole  trap  beam,  enough  to  selectively 
reflect  only  the  lattice  beam.  The  depth  of  the  lattice  is 
calibrated  using  Kapitza-Dirac  diffraction  of  a  23Na  BEC 
and  a  6Li2  molecular  BEC  and  by  lattice  modulation 
spectroscopy  on  the  6Li  cloud.  The  magnetic  field  and 
hyperfine  splitting  are  calibrated  using  rf  spectroscopy  on 
spin-polarized  clouds.  After  loading  the  lattice,  the  rf  pulse 
is  applied  for  a  duration  of  typically  1  ms.  Images  of  state 
1 2)  and  either  |1)  or  |3)  are  recorded  in  each  run  of  the 
experiment. 

To  ensure  loading  into  the  first  Bloch  band,  the  Fermi 
energy  and  temperature  of  the  cloud  are  kept  below  the 
second  band.  The  2D  Fermi  energy  E]n  =  l7I’'n" ,  with  n  the 
2D  density  per  spin  state,  is  typically  h  X  10  kHz.  The 
bottom  of  the  second  band  is  at  least  one  recoil  energy 
Er  =  h  X  29.3  kHz  above  the  bottom  of  the  first  band  in 
shallow  lattices  and  up  to  about  h  X  300  kHz  for  the  deep¬ 
est  lattices.  The  temperature  is  estimated  to  be  on  the  order 
of  the  Fermi  energy. 

rf  spectra  are  recorded  for  various  lattice  depths  and 
interaction  strengths.  Figure  1  shows  examples  of  spectra 
over  a  range  of  lattice  depths  at  the  3D  Feshbach  resonance 
and  on  the  BCS  side  of  the  resonance  at  721  G,  where 
fermion  pairing  in  3D  is  a  purely  many-body  effect.  At  the 
lowest  lattice  depths,  the  spectra  show  only  a  single  peak, 
shifted  to  positive  offset  frequencies  due  to  many-body 
interactions.  This  is  similar  to  the  case  without  a  lattice 
[24,30];  to  discern  a  peak  due  to  fermion  pairs  from  a  peak 
due  to  unbound  atoms  would  require  locally  resolved  rf 
spectroscopy  of  imbalanced  Fermi  gases  [30].  However,  as 
the  lattice  depth  is  raised,  the  single  peak  splits  into  two 
and  a  clear  pairing  gap  emerges.  The  narrow  peak  at  zero 
offset  is  the  free-to-free  transition,  and  the  asymmetric 
peak  at  positive  offset  is  the  pair  dissociation  spectrum. 
The  pair  spectrum,  especially  on  resonance,  shows  a  sharp 
threshold  and  a  long  tail  corresponding  to  dissociation  of 
fermion  pairs  into  free  atoms  with  nonzero  kinetic  energy. 

Binding  energies  are  determined  from  the  offset  fre¬ 
quency  of  the  pairing  threshold.  Although  the  line  shape 
in  Eq.  (2)  jumps  discontinuously  from  zero  to  its  maximum 
value,  the  spectra  are  observed  to  be  broadened.  This  is  to  a 
large  part  due  to  the  logarithmic  corrections  [28]  noted 
above,  which  predict  a  gradual  rise  at  the  threshold 
hv  =  Eb,  and  a  spectral  peak  that  is  slightly  shifted  from 
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RF  offset  (kHz) 

FIG.  1  (color  online).  Evolution  of  fermion  pairing  in  the 
3D-to-2D  crossover  in  a  one-dimensional  optical  lattice,  ob¬ 
served  via  rf  spectroscopy.  Shown  is  the  transferred  atom  num¬ 
ber  versus  rf  offset  frequency  relative  to  the  atomic  hyperfine 
splitting,  (a)  Spectra  at  the  Feshbach  resonance  at  690.7(1)  G 
with  d/ a  =  —0.01(4).  Lattice  depths  from  top  to  bottom  in  units 
of  Er :  1.84(3),  4.8(2),  6.1(2),  9.9(4),  12.2(4),  18.6(7),  and  19.5 
(7).  (b)  Spectra  on  the  BCS  side  at  720.7(1)  G,  d/ a  =  —  1.15(2). 
Lattice  depths  in  units  of  ER:  2.75(5),  4.13(7),  4.8(1),  6.0(2),  10.3 
(2),  and  18.1(4). 

Eh.  We  include  possible  additional  broadening  by  convolv¬ 
ing  the  theoretical  line  shape,  including  the  logarithmic 
correction,  with  a  Gaussian  function  of  width  wm.  The 
parameters  Eb  and  wm  are  determined  by  a  least-squares 
fit  to  the  measured  spectrum.  Typical  spectra  have  wm  of 
5  kHz,  consistent  with  our  estimates  of  broadening  based 
on  collisions  and  three-body  losses.  The  Fourier  broad¬ 
ening  is  1  kHz.  Power  broadening  is  about  5  kHz  on  the 
free-to-free  transition  and  less  than  1  kHz  on  the  bound-to- 
free  transition  due  to  the  reduced  wave  function  overlap. 
Inclusion  of  the  logarithmic  correction  is  found  to  be 
necessary  in  order  for  the  fit  function  to  reproduce  the 
observed  behavior  of  the  high-frequency  tail.  The  final- 
state  binding  energy  used  in  the  logarithmic  correction  for 
fitting  is  obtained  from  spectra  where  both  a  bound-to- 
bound  and  a  bound-to-free  peak  were  measured.  At  low 
lattice  depths,  the  2D  form  for  the  paired  spectrum  should 
differ  from  the  exact  shape  that  interpolates  between  the 
3D  and  2D  limits.  In  the  case  where  the  shape  of  the 
spectrum  is  given  by  the  3D  limit,  fitting  to  the  2D  form 
overestimates  the  binding  energy  by  8%. 

Figure  2  shows  the  measured  binding  energies  as  func¬ 
tion  of  V0/Er  for  several  interaction  strengths.  The  binding 


FIG.  2  (color  online).  Binding  energy  Eh  versus  lattice  depth 
V0  at  several  values  of  the  3D  scattering  length  a.  Eb  is 
normalized  via  the  lattice  frequency  coz.  Red  circles:  results 
from  spectra  at  690.7(1)  G  and  d/a  =  —0.01(4).  Green 
triangles:  720.7(1)  G,  d/a  =  —1.15(2).  Blue  squares: 
800.1(1)  G,  d/a  =  —2.69(1).  Curves  show  predictions  from 
Orso  et  al.  [11].  Horizontal  black  dashed  line:  harmonic  ap¬ 
proximation  result  for  l/a  =  0. 

energies  are  normalized  by  hwz  =  2-JV{]Er,  which  equals 
the  level  spacing  in  the  harmonic  approximation  to  the 
lattice  potential.  The  measured  binding  energies  grow 
with  increasing  lattice  depth  and  agree  reasonably  well 
with  theoretical  predictions  for  two-body  bound  pairs  in 
a  ID  lattice  [11].  The  binding  energy  at  the  3D  resonance 
approaches  a  constant  multiple  of  hwz  as  the  lattice 
depth  increases,  as  expected  from  the  2D  limit  [22,23]. 
Figure  3(a)  compares  the  binding  energies  measured  in 
lattices  deeper  than  11ER  to  predictions  in  the  harmonic 
quasi-2D  limit  [22,23].  At  the  3D  Feshbach  resonance,  we 
find  Eh  =  0.232(  1 6)hcoz  for  deep  lattices.  The  error  bar 
refers  to  the  standard  error  on  the  mean.  This  value  is  close 
to  the  harmonic  confinement  result  of  0. 244ft  w,  [23].  The 
exact  calculation  [11]  predicts  a  constant  downward  shift 
of  the  binding  energy  by  0.2 ER  for  deep  lattices  due  to  the 
anharmonicity  of  the  sinusoidal  potential.  For  V()  of  about 
20 Er,  this  gives  a  prediction  of  0.22 fuoz,  also  close  to  the 
measured  value. 

Figure  3(b)  shows  the  binding  energy  measured  in 
deep  lattices  normalized  by  the  exact  two-body  result 
[11]  versus  the  many-body  interaction  parameter 
\n{kFa2u)-  Overall,  the  binding  energies  are  close  to  the 
two-body  value,  even  in  the  strong-coupling  regime  for 
|  \n(kFa2o)\  <  1,  as  predicted  by  zero-temperature  mean- 
field  theory  [20].  The  data  show  a  slight  downward 
deviation  for  the  strongest  coupling.  At  fixed  reduced 
temperature  T/TF,  the  relationship  should  be  universal.  It 
will  thus  be  interesting  to  see  in  future  work  whether  the 
binding  energy  depends  significantly  on  temperature. 

The  bound-to-bound  transition  is  seen  in  Fig.  4  as  a 
naiTow  peak  at  negative  offset  frequencies.  In  the  regime 
where  Eb  can  be  found  from  the  pair  dissociation  spectrum. 
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FIG.  3  (color  online),  (a)  Binding  energy  of  fermion  pairs 
versus  interaction  strength  lj a  for  deep  lattices  (T0  >  17 ER). 
Solid  curve:  theoretical  prediction  in  the  2D  harmonic  limit 
[22,23].  (b)  Ratio  of  the  measured  binding  energy  to  the  two- 
body  result  [11]  versus  ln(kFa2r,)  for  V0>  \1ER.  Black  dia¬ 
monds:  binding  energy  determined  from  the  bound- to-bound 
transition  with  resonant  final-state  interactions.  Other  data  sym¬ 
bols:  see  Fig.  2.  Horizontal  line:  zero-temperature  mean-field 
theory  [20]. 


the  bound-to-bound  peak  position  directly  yields  the  bind¬ 
ing  energy  in  the  final  state  E'h.  For  example,  the  spectrum 
in  Fig.  4(a),  taken  at  the  3D  |1)  —  |3)  resonance  at 
690.7(1)  G  and  V0/ER  =  9.59(7),  gives  E'b/ER  =  18.0(1) 
at  a  final-state  interaction  of  d/a1  =  8.41(2).  Likewise,  the 
spectrum  in  Fig.  4(b)  at  Vq/Er  =  26.1(4)  and  a  magnetic 
field  of  751.1(1)  G,  where  d/a!  =  2.55(1),  gives  E'b/ER  = 
5.3(1).  An  independent  measurement  for  d/a  =  2.55(2) 
using  the  bound-to-free  spectrum  at  653.55  G  yields 
Eb/ER  =  5.25(2),  showing  that  bound-to-bound  transi¬ 
tions  correctly  indicate  binding  energies. 

The  BCS  side  of  the  2D  BEC-BCS  crossover  is  reached 
in  Fig.  4(c)  by  increasing  the  number  of  atoms  to  increase 
EF  and  increasing  the  magnetic  field  to  reach  a  lower 
binding  energy.  In  Fig.  4(c),  the  central  Fermi  energy  is 
h  X  43(6)  kHz  and  T/TF  =  0.5(2).  The  magnetic  field  is 
set  to  834.4(1)  G,  where  d/a  =  —3.06(1),  and  the  final- 
state  interactions  between  |1)  and  |2)  are  resonant,  with 
d/a'  =  —0.01(3).  The  lattice  depth  is  V0/ER  =  26.4(3). 
Thus,  we  know  that  E'b  =  0.232(16 )tuo,  =  2.4(2 )ER  at 
this  lattice  depth.  From  the  bound-to-bound  transition  in 
Fig.  4(c),  we  can  then  directly  determine  the  binding 
energy  of  |1)  —  |3)  fermion  pairs  to  be  Eh/ER  =  0.9(2). 
The  theoretical  prediction  [11]  for  two-body  binding  gives 
Eb/ER  =  0.82(1).  The  measured  binding  energy  gives  a 
many-body  interaction  parameter  of  ln(kfa2o)  =  0.6(1), 
on  the  BCS  side  but  within  the  strongly  interacting  regime, 
where  one  expects  many-body  effects  beyond  mean-field 
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FIG.  4  (color  online).  Spectra  including  the  bound- to-bound 
transition,  a  narrow  peak  at  negative  rf  offset.  Shown  are  spectra 
at  magnetic  fields  of  (a)  690.7(1)  G,  (b)  751.1(1)  G,  and 
(c)  834.4(1).  The  interaction  parameters  d/a  are  (a)  —0.01(4), 
(b)  —1.91(1),  and  (c)  —3.06(1).  Lattice  depths  in  units  of  ER  are 
(a)  9.59(7),  (b)  26.1(4),  and  (c)  26.4(3).  The  bound-to-free 
transition  is  not  visible  in  (c).  The  transfer  is  from  |1)  to  |2) 
in  (a)  and  (b)  and  from  |3)  to  |2)  in  (c). 


BEC-BCS  theory  [26,29].  It  is  therefore  interesting  that  the 
measured  binding  energy  is  close  to  the  expected  two-body 
binding  energy  to  much  better  than  the  Fermi  energy,  as 
predicted  by  mean-field  theory  [20]. 

In  conclusion,  we  have  measured  the  binding  energy  of 
fermion  pairs  along  the  crossover  from  3D  to  2D  in  a  one¬ 
dimensional  optical  lattice.  Measurements  were  performed 
at  several  lattice  depths  and  scattering  lengths,  allowing 
quantitative  comparison  with  theoretical  predictions. 
Considering  the  fact  that  the  gas  is  a  strongly  interacting 
many-body  system,  the  close  agreement  with  two-body 
theory  is  surprising,  especially  in  the  strong-coupling  re¬ 
gime.  While  mean-field  BEC-BCS  theory  in  2D  predicts 
this  behavior  [20],  it  misses  other  important  features  of  the 
many-body  system,  most  strikingly  the  interaction  between 
fermion  pairs  [13].  Superfluidity  in  a  one-dimensional 
lattice  will  be  an  exciting  topic  for  future  studies.  Stacks 
of  weakly  coupled,  superfluid  2D  layers  would  constitute  a 
basic  model  of  the  geometry  found  in  high-temperature 
superconductors. 
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Feynman  diagrams  versus  Fermi-gas 
Feynman  emulator 

K.  Van  Houcke1,2*,  F.  Werner1,3,  E.  Kozik4,5,  N.  Prokof'ev1,6,  B.  Svistunov1,6,  M.  J.  H.  Ku7, 
A.  T.  Sommer7,  L.  W.  Cheuk7,  A.  Schirotzek8  and  M.  W.  Zwierlein7 


Precise  understanding  of  strongly  interacting  fermions,  from 
electrons  in  modern  materials  to  nuclear  matter,  presents 
a  major  goal  in  modern  physics.  However,  the  theoretical 
description  of  interacting  Fermi  systems  is  usually  plagued 
by  the  intricate  quantum  statistics  at  play.  Here  we  present 
a  cross-validation  between  a  new  theoretical  approach,  bold 
diagrammatic  Monte  Carlo1-3,  and  precision  experiments  on 
ultracold  atoms.  Specifically,  we  compute  and  measure,  with 
unprecedented  precision,  the  normal-state  equation  of  state  of 
the  unitary  gas,  a  prototypical  example  of  a  strongly  correlated 
fermionic  system4-6.  Excellent  agreement  demonstrates  that  a 
series  of  Feynman  diagrams  can  be  controllably  resummed  in  a 
non-perturbative  regime  using  bold  diagrammatic  Monte  Carlo. 

In  his  seminal  1981  lecture7,  Richard  Feynman  argued  that 
an  arbitrary  quantum  system  cannot  be  efficiently  simulated 
with  a  classical  universal  computer,  because  generally,  quantum 
statistics  can  only  be  imitated  with  a  classical  theory  if  probabilities 
are  replaced  with  negative  (or  complex)  weighting  factors.  For 
the  majority  of  many-particle  models  this  indeed  leads  to  the 
so-called  sign  problem,  which  has  remained  an  insurmountable 
obstacle.  According  to  Feynman,  the  only  way  out  is  to  employ 
computers  made  out  of  quantum-mechanical  elements7.  The  recent 
experimental  breakthroughs  in  cooling,  probing  and  controlling 
strongly  interacting  quantum  gases  prompted  a  challenging 
effort  to  use  this  new  form  of  quantum  matter  to  realize 
Feynman’s  emulators  of  fundamental  microscopic  models7,8. 
Somewhat  ironically,  Feynman’s  arguments,  which  led  him  to  the 
idea  of  emulators,  may  be  defied  by  a  theoretical  method  that 
he  himself  devised,  namely  Feynman  diagrams.  This  technique 
organizes  the  calculation  of  a  given  physical  quantity  as  a  series  of 
diagrams  representing  all  the  possible  ways  particles  can  propagate 
and  interact  (for  example,  ref.  9).  For  the  many-body  problem,  this 
diagrammatic  expansion  is  commonly  used  either  in  perturbative 
regimes  or  within  uncontrolled  approximations.  However,  the 
introduction  of  diagrammatic  Monte  Carlo  recently  allowed  one  to 
go  well  beyond  the  first  few  diagrams,  and  even  reach  convergence 
of  the  series  in  a  moderately  correlated  regime1,10. 

In  this  Letter  we  show  that  for  a  strongly  correlated  system  and 
down  to  a  phase  transition,  the  diagrammatic  series  can  still  be 
given  a  mathematical  meaning  and  leads  to  controllable  results 
within  bold  diagrammatic  Monte  Carlo  (BDMC).  This  approach, 
proposed  in  refs  1-3,  is  first  implemented  here  for  the  many-body 
problem.  We  focus  on  the  unitary  gas,  that  is,  spin- 1/2  fermions 


with  zero-range  interactions  at  infinite  scattering  length4-6.  This 
system  offers  the  unique  possibility  to  stringently  test  our  theory 
against  a  quantum  emulator  realized  here  with  trapped  ultracold 
6Li  atoms  at  a  broad  Feshbach  resonance4-6.  This  experimental 
validation  is  indispensable  for  our  theory,  based  on  resummation  of 
a  possibly  divergent  series:  although  the  physical  answer  is  shown  to 
be  independent  of  the  applied  resummation  technique — suggesting 
that  the  procedure  is  adequate — its  mathematical  validity  remains 
to  be  proven.  In  essence,  nature  provides  the  ‘proof.  This  presents 
the  first — although  long-anticipated — compelling  example  of  how 
ultracold  atoms  can  guide  new  microscopic  theories  for  strongly 
interacting  quantum  matter. 

At  unitarity,  the  disappearance  of  an  interaction-imposed  length 
scale  leads  to  scale  invariance.  This  property  renders  the  model 
relevant  for  other  physical  systems  such  as  neutron  matter.  It 
also  makes  the  balanced  (that  is,  spin- unpolarized)  unitary  gas 
ideally  suited  for  the  experimental  high-precision  determination 
of  the  equation  of  state  (EOS)  described  below.  Finally,  it 
implies  the  absence  of  a  small  parameter,  making  the  problem 
notoriously  difficult  to  solve. 

In  traditional  Monte  Carlo  approaches,  which  simulate  a 
finite  piece  of  matter,  the  sign  problem  causes  an  exponential 
increase  of  the  computing  time  with  system  size  and  inverse 
temperature.  In  contrast,  BDMC  simulates  a  mathematical  answer 
in  the  thermodynamic  limit.  This  radically  changes  the  role  of  the 
fermionic  sign.  Diagrammatic  contributions  are  sign-alternating 
with  order,  topology  and  values  of  internal  variables.  Because 
the  number  of  graphs  grows  factorially  with  diagram  order,  a 
near-cancellation  between  these  contributions  is  actually  necessary 
for  the  series  to  be  resummable  by  techniques  requiring  a 
finite  radius  of  convergence.  We  find  that  this  ‘sign  blessing’ 
indeed  takes  place. 

In  essence,  BDMC  solves  the  full  quantum  many-body  problem 
by  stochastically  summing  all  the  skeleton  diagrams  for  irreducible 
single-particle  self-energy  E  and  pair  self-energy  77,  expressed 
in  terms  of  bold  (that  is,  fully  dressed)  single-particle  and  pair 
propagators  G  and  F  which  are  determined  self-consistently 
(Fig.  1 ) .  The  density  EOS  (that  is,  the  relation  between  total  density 
n,  chemical  potential  /z  and  temperature  T)  is  given  by  G  at 
zero  distance  and  imaginary  time,  n(/z,  T)  =  2  G(r  =  0,  r  =  0-). 
The  thermodynamic  limit  can  be  taken  analytically.  The  sum  of 
ladder  diagrams  built  on  the  bare  single-particle  propagator  defines 
a  partially  dressed  pair  propagator  F°.  As  F°  is  well  defined 
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Figure  1  |  Bold  diagrammatic  Monte  Carlo  The  skeleton  diagrammatic 
series  for  the  self-energy  E  and  the  pair  self-energy  77  is  evaluated 
stochastically  (lower  box).  The  diagrams  are  built  on  dressed  one-body 
propagators  6  and  pair  propagators  r,  which  themselves  are  the  solution 
of  the  Dyson  and  Bethe-Salpeter  equations  (upper  box).  This  cycle  is 
repeated  until  convergence  is  reached.  G°  is  the  non-interacting  propagator 
and  r°  is  the  partially  dressed  pair  propagator  obtained  by  summing  the 
bare  ladder  diagrams. 

for  the  zero-range  continuous-space  interaction,  the  zero-range 
limit  can  also  be  taken  analytically.  This  is  in  sharp  contrast  with 
other  numerical  methods11-13,  where  taking  the  thermodynamic 
and  zero-range  limits  is  computationally  very  expensive.  BDMC 
performs  a  random  walk  in  the  space  of  irreducible  diagrams  using 
local  updates.  The  simulation  is  run  in  a  self-consistent  cycle  (along 
the  lines  of  ref.  2)  until  convergence  is  reached.  Full  details  will 
be  presented  elsewhere.  In  essence,  our  approach  upgrades  the 
standard  many-body  theories  based  on  one  lowest-order  diagram 
(for  example,  refs  14,15)  to  millions  of  graphs. 

In  the  quantum  degenerate  regime,  we  do  not  observe 
convergence  of  the  diagrammatic  series  for  E  and  77  evaluated 
up  to  order  9.  Here,  order  N  means  17-diagrams  with  N  vertices 
(that  is,  N  F-lines)  and  77-diagrams  with  N  —  1  vertices.  To 
extract  the  infinite-order  result,  we  apply  the  following  Abelian 
resummation  methods16.  The  contribution  of  all  diagrams  of  order 
N  is  multiplied  by  e~eXN-' ,  where  A„  depends  on  the  resummation 
method:  (1)  A„  =  n  logn  (with  A0  =  0)  for  Lindelof16,  (2)  A„  =  (n—  1) 
log(n—  1)  (with  A0  =  Ar  =  0)  for  ‘shifted  Lindelof,  or  (3)  A„  =  n2  for 
Gaussian17.  A  full  simulation  is  performed  for  each  e,  and  the  final 
result  is  obtained  by  extrapolating  to  e  =  0  (Fig.  2). 

This  protocol  relies  on  the  following  crucial  mathematical 
assumptions:  (1)  the  JVth  order  contribution  of  the  diagrammatic 
expansion  for  E  (for  fixed  external  variables)  is  the  Nth  coefficient 
of  the  Taylor  series  at  z  =  0  of  a  function  g(z)  which  has  a  non-zero 
convergence  radius,  (2)  the  analytic  continuation  ^(1),  performed 
by  the  above  resummation  methods16,17,  is  the  physically  correct 
value  of  E.  The  same  assumptions  should  hold  for  77. 

Proving  these  assumptions  is  an  open  mathematical  challenge. 
Note  that  Dyson’s  collapse  argument18  is  not  applicable  to  immedi¬ 
ately  disprove  the  assumption  ( 1 )  of  a  non-zero  convergence  radius: 
indeed,  unlike  QED,  our  skeleton  series  is  not  an  expansion  in 
powers  of  a  coupling  constant  whose  sign  change  would  lead  to 
an  instability.  The  first  important  evidence  for  the  validity  of  our 
mathematical  assumptions  is  that  the  three  different  resummation 
methods  yield  consistent  results.  For  an  independent  test,  we 
turn  to  experiments. 

The  present  experiment  furnishes  high-precision  data  for  the 
density  n  as  a  function  of  the  local  value  V  of  the  trapping  potential 
(Fig.  3  and  Methods).  We  start  the  process  by  obtaining  the  EOS  at 
high  temperatures  in  the  non-degenerate  wings  of  the  atom  cloud, 


Figure  2  |  Cross-validation  between  resummation  procedure  and 
experiment  at  =  +1.  Bold  diagrammatic  Monte  Carlo  data  for  the 
dimensionless  density  nA3,  as  a  function  of  the  parameter  e  controlling  the 
resummation  procedure,  for  three  different  resummation  methods:  Lindelof 
(blue  circles),  shifted  Lindelof  (black  diamonds),  and  Gauss  (open  green 
squares).  The  solid  lines  are  linear  fits  to  the  Monte  Carlo  data,  their  e-»0 
extrapolation  agrees  within  error  bars  with  the  experimental  data  point 
(filled  red  square).  (In  the  opposite  limit  e  — x  oo,  the  Lindelof  (resp.  shifted 
Lindelof)  curves  will  asymptote  to  the  first15,21  (resp.  third)  order  results, 
shown  by  the  dashed  (resp.  dash-dotted)  line.)  Error  bars  for  each  e 
represent  the  statistical  error,  together  with  the  estimated  systematic  error 
coming  from  not  sampling  diagrams  of  order  >9. 

where  the  virial  expansion  is  applicable.  Once  the  temperature 
and  the  chemical  potential  have  been  determined  from  fits  to  the 
wings  of  the  cloud,  the  data  closer  to  the  cloud  centre  provides 
a  new  prediction  of  the  EOS.  The  process  is  iterated  to  access 
lower  temperatures. 

Scale  invariance  allows  one  to  write  the  density  EOS  as 
n(/z,  T)  A3  =/ (/3/rx),  with  A  =  V2nh2/{mkBT)  the  thermal  de  Broglie 
wavelength,  ft  =  l/(fcBT)  the  inverse  temperature  and /  a  universal 
function.  A  convenient  normalization  of  the  data  is  provided 
by  the  EOS  of  a  non-interacting  Fermi  gas,  n0 A3  =f0(ftii).  In 
Fig.  4a,  we  thus  report  the  ratio  «(/z,  T)/n0(fJL,  T )  =/(y6/z)//o(JS/z), 
bringing  out  the  difference  between  the  ideal  and  the  strongly 
interacting  Fermi  gas.  The  Gibbs-Duhem  relation  allows  us  to  also 
calculate  the  pressure  at  a  given  chemical  potential,  P(/z0,  T)  = 
fZe d/x  n(/i,  T)  =  \/(ft\2)F(ftnQ),  where  F{x)  =  f^dx’f  (xr).  We 
normalize  it  by  the  pressure  of  the  ideal  Fermi  gas  and  show 
F(ft[i)/F0(ft[i)  (Fig.  4b).  The  agreement  between  BDMC  and 
experiment  is  excellent.  The  comparison  is  sufficiently  sensitive  to 
validate  the  procedure  of  resumming  and  extrapolating  (Fig.  2). 
The  result  was  checked  to  be  independent  of  the  maximal  sampled 
diagram  order  Nmax  e  {7;  8;  9)  within  the  error  bars  shown  in 
Fig.  2  for  each  e.  The  BDMC  final  error  bar  in  Fig.  4  is  the 
sum  of  the  conservatively  estimated  systematic  errors  from  the 
uncertainty  of  the  e  — >  0  extrapolation  and  from  the  dependence 
on  numerical  grids  and  cutoffs,  the  latter  being  reduced  by 
analytically  treating  high-momentum  short-time  singular  parts. 
The  systematic  error  in  the  experiment  is  determined  to  be 
about  1%  by  the  independent  determination  of  the  EOS  of 
the  non-interacting  Fermi  gas.  The  experimental  error  bars  of 
Fig.  4  also  include  the  statistical  error,  which  is  <0.5%,  thanks 
to  the  scale  invariance  of  the  balanced  unitary  gas:  irrespective 
of  shot-to-shot  fluctuations  of  atom  number  and  temperature, 
all  experimental  profiles  contribute  to  the  same  scaled  EOS- 
function  /.  The  dominant  uncertainty  on  the  experimental  EOS 
stems  from  the  uncertainty  in  the  position  of  the  6Li  Feshbach 
resonance,  known  to  be  at  834.15  ±  1.5  G  from  spectroscopic 
measurements19.  The  change  in  energy,  pressure  and  density  with 
respect  to  the  interaction  strength  is  controlled  by  the  so-called 
contact20  that  is  obtained  from  F  in  the  BDMC  calculation. 
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Figure  3  |  Constructing  the  EOS  from  in  situ  imaging.  The  atom  cloud  shown  contains  N  =  8  x  104  atoms  for  each  spin  state,  with  a  local  Fermi  energy  of 
Ef  =  370  nK  at  the  centre,  a,  Absorption  image  of  the  atomic  cloud  after  quadrant  averaging,  b,  Reconstructed  local  density  n(p,z).  c,  Equipotential 
averaging  produces  a  low-noise  density  profile,  n  versus  V.  Thermometry  is  performed  by  fitting  the  experimental  data  (red)  to  the  known  portion  of  the 
EOS  (solid  blue  line),  starting  with  the  virial  expansion  for  <  —1.25  (green  dashed  line).  In  this  example,  the  EOS  is  known  for  ft/j.  <  —  0.25,  and  the  fit  to 
the  density  profile  yields  T  =  113  nK,  and  /i/c  =  1.63.  d,  Given  p  and  T,  the  density  profile  can  be  rescaled  to  produce  the  EOS  nA3  versus  /J/z. 


p/(cBT  ju/kBT 


Figure  4  |  Equation  of  state  of  the  unitary  Fermi  gas  in  the  normal  phase.  Density  n  (a)  and  pressure  P  (b)  of  a  unitary  Fermi  gas,  normalized  by  the 
density  no  and  the  pressure  Pq  of  a  non-interacting  Fermi  gas,  versus  the  ratio  of  chemical  potential  /z  to  temperature  T.  Blue  filled  squares:  BDMC  (this 
work),  red  filled  circles:  experiment  (this  work).  The  BDMC  error  bars  are  estimated  upper  bounds  on  systematic  errors.  The  error  bars  are  one  standard 
deviation  systematic  plus  statistical  errors,  with  the  additional  uncertainty  from  the  Feshbach  resonance  position  shown  by  the  upper  and  lower  margins  as 
red  solid  lines.  Black  dashed  line  and  red  triangles:  Theory  and  experiment  (this  work)  for  the  ideal  Fermi  gas,  used  to  assess  the  experimental  systematic 
error.  Green  solid  line:  third  order  virial  expansion.  Open  squares:  first  order  bold  diagram15,21.  Green  open  circles:  Auxiliary  Field  QMC  (ref.  11).  Star: 
superfluid  transition  point  from  Determinental  Diagrammatic  Monte  Carlo13.  Filled  diamonds:  experimental  pressure  EOS  (ref.  22).  Open  pentagons: 
pressure  EOS  (ref.  23). 
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This  allows  us  to  define  the  uncertainty  margins  above  and  below 
the  experimental  data  (Fig.  4)  that  give  the  prediction  for  the  unitary 
EOS  if  the  true  Feshbach  resonance  lies  1.5  G  below  or  above 
834.15  G,  respectively. 

Our  results  clearly  differ  from  previous  theoretical  and  ex¬ 
perimental  results.  Deviations  from  the  theory  based  on  the 
first-order  Feynman  diagrams15,21  are  expected,  and  rather  re¬ 
markably  moderate.  Differences  with  lattice  Monte  Carlo  data11,13 
may  seem  more  surprising,  as  in  the  particular  case  of  the  bal¬ 
anced  system  these  algorithms  are  free  of  the  sign  problem,  al¬ 
lowing  one  in  principle  to  approach  the  balanced  unitary  gas 
model  in  an  unbiased  way.  However,  eliminating  systematic  er¬ 
rors  from  lattice-discretization  and  finite  volume  requires  extrap¬ 
olations  which  are  either  not  done11  or  difficult  to  control12,13. 
The  ENS  experimental  pressure  EOS  (ref.  22)  lies  systemati¬ 
cally  below  ours,  slightly  outside  the  reported  error  bar.  The 
experimental  results  from  Tokyo23  do  not  agree  with  the  virial 
expansion  at  high  temperature.  The  BDMC  results  agree  well 
with  the  present  experimental  data  all  the  way  down  to  the 
critical  temperature  for  superfluidity  (Fig.  4).  On  approaching 
(Pi u)0  we  observe  the  growth  of  the  correlation  length  in  the 
BDMC  pair  correlation  function  r.  A  protocol  for  extracting 
the  critical  temperature  itself  from  the  BDMC  simulation  will 
be  presented  elsewhere. 

We  are  not  aware  of  any  system  of  strongly  correlated  fermions 
in  nature  where  experimental  and  unbiased  theoretical  results 
were  compared  at  the  same  level  of  accuracy.  Even  for  bosons, 
the  only  analogue  is  liquid  4He.  This  promotes  the  unitary  gas 
to  the  major  testing  ground  for  unbiased  theoretical  treatments. 
The  present  BDMC  implementation  should  remain  applicable  at 
finite  polarization  and/or  finite  scattering  length,  opening  the 
way  to  rich  physics  which  was  already  addressed  by  cold  atoms 
experiments6,24-27.  We  also  plan  to  extend  BDMC  to  superfluid 
phases  by  introducing  anomalous  propagators.  Moreover,  as 
the  method  is  generic,  we  expect  numerous  other  important 
applications  to  long-standing  problems  across  many  fields. 

Note  added  in  proof:  After  a  preprint  of  this  work  became  available, 
new  auxiliary-field  quantum  Monte  Carlo  data  were  presented28, 
with  undetermined  systematic  errors  whose  evaluation  in  future 
work  is  called  for  by  the  authors  of  ref.  28. 

Methods 

The  experimental  set-up  has  been  described  previously24.  In  short,  ultracold 
fermionic  6  Li  is  brought  to  degeneracy  by  sympathetic  cooling  with  23Na.  A 
two-state  mixture  of  the  two  lowest  hyperfine  states  of  6Li  is  further  cooled  in  a 
hybrid  magnetic  and  optical  trap  at  the  broad  Feshbach  resonance  at  834  G.  We 
employ  high-resolution  in  situ  absorption  imaging  to  obtain  the  column  density  of 
the  gas,  that  is  converted  into  the  full  3D  density  using  the  inverse  Abel  transform29. 
Equidensity  lines  provide  equipotential  lines  that  are  precisely  calibrated  using 
the  known  axial,  harmonic  potential  (axial  frequency  vz  =  22.83  ± 0.05  Hz). 
Equipotential  averaging  yields  low-noise  profiles  of  density  n  versus  potential  V. 
Density  is  absolutely  calibrated  by  imaging  a  highly  degenerate,  highly  imbalanced 
Fermi  mixture,  and  fitting  the  majority  density  profile  to  the  ideal  Fermi  gas  EOS 
(ref.  24).  In  contrast  to  previous  studies22,23,  our  analysis  does  not  rely  on  the 
assumption  of  harmonic  trapping. 

Thermometry  is  performed  by  fitting  the  density  profile  to  the  EOS 
constructed  thus  far,  restricting  the  fit  to  the  portion  of  the  density  profile 
where  the  EOS  is  valid.  In  the  high-temperature  regime,  the  EOS  is  given  by 
the  virial  expansion 

nX3  =  +  2  b2e2^  +  3  b3e3^  +  •  •  ■  (1) 

where  the  virial  coefficients  are  b2  =  3^2/8  (ref.  30),  and  b3  =  —0.29095295 
(ref.  31).  Fitting  a  high-temperature  cloud  to  the  virial  expansion  gives  the 
temperature  T  and  the  chemical  potential  [i  of  the  cloud,  and  the  EOS  nX3  =f  {fill) 
can  be  constructed.  We  have  used  equation  (1)  for  /J/z  <  (fill) max  =  — 1.25  and  we 
checked  that  our  EOS  did  not  change  within  statistical  noise  if  we  instead  used 
(/J/z)max  =  —0.85.  Once  a  new  patch  of  EOS  has  been  produced,  it  can  then  in  turn 
be  used  to  fit  colder  clouds.  Iteration  of  this  method  allows  us  to  construct  the  EOS 
to  arbitrarily  low  temperature.  A  total  of  ~  1,000  profiles  were  used,  with  10-100 
profiles  (50  on  average)  contributing  at  any  given  f/i. 
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